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SUMMARY 

A model  has  heen  formulated  which,  when  completed,  will  represent 
the  main  components  and  interactions  that  are  believed  to  be  important 
in  the  ecosystem  of  Lake  Conway  and  other  southeastern  lakes  as  well. 

The  model  includes  phytoplankton,  epipelic  algae  (i.e.  that  living  upon 
the  sediment),  and  macrophytic  plants  with  their  associated  epiphytic 
algae,  as  well  as  zooplankton,  benthic  invertebrates,  herbivorous  fish, 
young  and  adult  primary  predator  fish,  and  secondary  predator  fish.  In 
the  model,  the  lake  is  divided  into  epilimnion  and  hypolimnion,  and  the 
sediments  comprise  a third  physical  unit. 

Information  from  both  field  work  and  literature  surveys  was  used 

in  formulating  the  model.  Productivity  studies  on  the  lake  provided  an 

2 

estimate  of  gross  primary  productivity  of  898  g C/m  per  year,  and 
community  respiration  was  nearly  equal  to  it.  Changes  in  biomass  of 
components  were  used  that  were  reported  by  other  researchers . Tempera- 
ture measurements  made  during  the  study  year  were  used  in  modeling  the 
feeding  and  respiration  functions  of  various  components.  Shading  effects 
of  the  phytoplankton  and  macrophytic  plants  on  each  other  and  on 
epipelic  algae  were  taken  into  account  in  programming  seasonal  changes 
in  photosynthesis.  Macrophytic  plants  were  assumed  to  obtain  nutrients 
from  the  sediments  and  to  release  nutrients  into  the  surrounding  water 
at  a rate  proportional  to  their  rate  of  respiration.  The  model  as 
currently  formulated  requires  better  definition  of  nutrient-productivity 
relationships  and  of  zooplankton  population  dynamics. 

Given  simulated  introduction  of  white  amur,  the  model  predicts  de- 
creases in  phytoplankton,  zooplankton,  and  herbivorous  fish  and  in- 
creases in  epipelic  algae,  benthic  invertebrates,  and  secondary  predator 
fish,  all  for  a 10-year  period.  The  model  also  predicts  a decrease  in 
the  quantity  of  sediments,  which  suggests  a possible  reversal  of  trend 
in  the  trophic  state  of  the  lake.  However,  these  results  may  change 
with  further  refinement  of  the  model. 
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PREFACE 


The  work  described  in  this  volume  was  performed  under  Contract 
No.  DACW39-76-C-0019  between  the  U.  S.  Army  Engineer  Waterways  Experiment 
Station  (WES),  Vicksburg,  Miss.,  and  the  University  of  Florida, 
Gainesville.  The  work  was  sponsored  by  the  U.  S.  Army  Engineer  District, 
Jacksonville,  and  by  the  Office,  Chief  of  Engineers,  U.  S.  Army. 

This  is  the  seventh  of  eight  volumes  that  constitute  the  first  of 
a series  of  reports  documenting  a large-scale  operations  management  test 
of  use  of  the  white  amur  for  control  of  problem  aquatic  plants  in  Lake 
Conway,  Fla.  Report  1 presents  the  results  of  the  baseline  studies  of 
Lake  Conway;  subsequent  reports  will  present  the  annual  poststocking 
results. 

The  work  was  performed  and  this  volume  was  written  by 
Dr.  Katherine  C.  Ewel  and  Mr.  Thomas  D.  Fontaine  of  the  School  of 
Forest  Resources  and  Conservation  of  the  University  of  Florida. 

Mr.  Ronald  L.  Myers  and  Mr.  James  Sampson  assisted  at  several  stages  of 
the  project.  Mr.  Glenn  Smerage  and  Dr.  Jerome  Shireman  made  useful 
comments  on  an  earlier  version  of  the  model,  and  Dr.  H.  T.  Odom  made 
several  helpful  suggestions  on  model  development  and  productivity  data 
interpretation.  The  Shenandoah  Park  Residents'  Association  at  Lake 
Conway  generously  allowed  the  authors  use  of  their  dock  and  launching 
ramp  for  the  productivity  study. 

The  work  was  monitored  at  WES  in  the  Mobility  and  Environmental 
Systems  Laboratory  (MESL)  under  the  general  supervision  of  Mr.  W.  G. 
Shockley,  Chief  of  MESL,  and  Mr.  B.  0.  Benn,  Chief  of  the  Environmental 
Systems  Division,  and  under  the  direct  supervision  of  Mr.  J.  L.  Decell, 
Chief  of  the  Aquatic  Plant  Research  Branch  (APRB).  The  APRB  is  now  part 
of  the  recently  organized  Environmental  Laboratory  of  which  Dr.  John 
Harrison  is  Chief. 

Director  of  WES  during  the  period  of  the  contract  was  COL  J.  L. 
Cannon,  CE.  Technical  Director  was  Mr.  F.  R.  Brown. 
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LARGE-SCALE  OPERATIONS  MANAGEMENT  TEST  OF  USE  OF  THE  WHITE 
AMUR  FOR  CONTROL  OF  PROBLEM  AQUATIC  PLANTS 

P'-SELINE  STUDIES 

A Model  for  Evaluation  of  the  Response  of  the 
Lake  Conway,  Florida,  Ecosystem  to 
Introduction  of  the  White  Amur 


PART  I : INTRODUCTION 

Background 

1.  The  U.  S.  Army  Corps  of  Engineers  is  currently  conducting  a 
large-scale  field  test  of  the  ability  of  the  white  amur  ( Ctenopharyngodon 
i della)  to  control  excessive  growths  of  hydrilla  (Hydrilla  verticillata) 
in  Lake  Conway,  near  Orlando,  Fla.  Lake  Conway  comprises  a chain  of  five 
bodies  of  water  and  has  a surface  area  of  7-6  km  (Figure  l).  Maximum 
depth  varies  from  8 to  12  m among  the  five  pools,*  but  the  mean  depth  of 
all  pools  is  approximately  5m.  As  a result  of  the  lake's  lobed  mor- 
phology, it  has  2.6  times  the  shoreline  that  it  would  have  if  it  were 
circular  but  had  the  same  area.  Because  of  this,  it  has  considerable 
potential  for  littoral  zone  plant  proliferation.  Since  the  lakeshore 
has  been  developed  for  residential  use,  however,  only  small  areas  of 
emergent  macrophytes  still  remain.  Nevertheless,  submersed  plants  are 
common  in  depths  of  up  to  4 m.  Considerable  stands  of  pondweed 
(Potamogeton  illinoensis)  are  found  in  shallow  areas,  and  mats  of 
hydrilla  have  covered  large  areas  of  open  water  in  the  past.  During  the 
study  year,  1976-1977*  Nitella  megacarpa  was  the  main  submerged  aquatic 
plant  found  in  Lake  Conway. 

Scope 

2.  This  report  describes  the  development  and  summarizes  the 
* E.  Blancher,  personal  communications. 
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Figure  1.  Map  of  Lake  Conway,  Fla.,  showing  locations 
of  transects  used  for  productivity  measurements 

content  of  a mathematical  model  that  has  been  formulated  in  order  to 
predict  the  long-term  effects  of  introduction  of  the  white  amur  not  only 
on  the  population  of  hydrilla  but  on  the  entire  Lake  Conway  ecosystem  as 
well.  In  addition  to  the  modeling  effort,  community  and  planktonic 
metabolism  data  were  collected  on  Lake  Conway  for  1 year,  1976-1977- 
These  data  indicate  the  relative  importance  of  different  producer  popula- 
tions in  the  lake.  Data  on  levels  of  producer  and  consumer  populations 
as  well  as  water  chemistry,  which  were  collected  by  researchers  working 
on  other  segments  of  the  baseline  studies,  are  included  in  the  model. 


PART  II:  METABOLISM  MEASUREMENTS 


Methods 

Temperature  and  dissolved  oxygen 

3.  Measurements  of  temperature  and  dissolved  oxygen  were  made 
monthly  at  three  stations  on  each  of  three  transects  in  the  East  Pool  of 
Little  Lake  Conway  (Figure  l).  The  three  transects  are  each  about  100  m 
long  and  extend  from  shorelines  having  different  degrees  of  development. 
Transect  1 extends  from  a forested  area  but  runs  parallel  to  a resi- 
dential shoreline.  Transect  2,  located  near  the  joining  of  the  two 
pools  comprising  Little  Lake  Conway,  is  influenced  primarily  by  a for- 
ested shoreline  and  cattail  marsh.  Transect  3,  which  is  most  heavily 
influenced  by  urban  development,  is  across  the  lake  from  Transect  2.  It 
starts  at  a lawn  and  is  underlain  by  sand  for  part  of  its  length. 

1+.  The  three  stations  at  each  transect  were  characterized  by  dif- 
ferent vegetational  communities.  Pondweed  and  eelgrass  (Vallisneria 
americana)  were  found  at  the  most  shallow  stations  (2  m);  hydrilla  and 
nitella  were  characteristic  of  the  medium-depth  stations  (U  m);  and 
phytoplankton  were  dominant  at  the  pelagic,  deepwater  stations 
(6  to  8 m) . 

5.  Water  samples  for  dissolved  oxygen  determination  were  col- 
lected with  a Van  Dorn  bottle  and  analyzed  by  the  Winkler  method. 

Samples  were  taken  at  the  surface,  at  0.5  m,  from  1 to  6 m at  1-m  inter- 
vals at  the  two  deepest  stations,  and  at  0.5-m  intervals  at  the  marsh. 
Samples  were  collected  as  close  to  dawn  and  dusk  as  possible,  with  an 
additional  sample  in  the  middle  of  the  night  and  from  one  to  five  addi- 
tional daytime  samples. 

6.  Light  and  dark  bottles  were  also  incubated  each  month  in  the 
photic  zone  at  each  station  on  one  of  the  transects  in  order  to  differ- 
entiate between  planktonic  and  macrophytic  contributions  to  the  produc- 
tion and  respiration  rates  of  the  entire  community.  Incubation  time  was 
generally  from  U to  6 hours  during  the  day. 

7.  On  three  occasions,  under  conditions  ranging  from  smooth  to 
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choppy  water,  measurements  of  the  rate  of  diffusion  of  oxygen  from  the 

water  to  the  air  were  made.  A plastic  diffusion  dome  painted  gray  (to 

stop  internal  heat  buildup)  was  used.  The  internal  volume  of  the  dome 

2 

was  5«T  £ and  the  surface  area  was  0.1225  m . The  dome  was  purged  with 
nitrogen  gas  until  an  oxygen  probe  inside  the  dome  (previously  cali- 
brated to  100  percent  saturation  in  the  air)  registered  zero  percent 
saturation.  The  rate  of  change  of  percent  oxygen  saturation  was  then 
monitored. 

8.  In  addition  to  the  three  transects  described  above,  two  other 
transects  were  established — one  in  each  of  the  other  two  large  pools  of 
Lake  Conway  (Figure  l).  Dissolved  oxygen  was  measured  with  an  oxygen 
meter  at  stations  along  these  transects.  Results  after  several  months 
indicated  that  the  differences  between  these  stations  and  the  ones  in  the 
East  Pool  were  not  significant  enough  to  warrant  continuing  these  addi- 
tional measurements. 

Calculations  of  diurnal  curves 

9.  To  calculate  productivity  on  an  areal  basis  from  the  dis- 
solved oxygen  measurements,  it  was  necessary  to  assess  dissolved  oxygen 
changes  in  1-m-deep  layers  of  water.  Dissolved  oxygen  levels  at  the 
surface,  at  0.5  m,  and  at  1 m were  averaged  to  obtain  a representative 
value  for  the  first  metre-deep  layer.  Measurements  at  1 and  2 m were 
then  averaged  to  obtain  a value  for  the  second  layer,  and  so  on. 

10.  An  example  of  a diurnal  oxygen  curve  is  given  for  Station  IB 
in  Figure  2.  In  this  case,  the  oxygen  concentration  in  the  water  column 
reached  its  highest  value  approximately  2 hours  before  sundown  and  then 
declined  to  its  lowest  value  at  dawn.  Surface  temperature  and  percent 
oxygen  saturation  are  also  plotted  so  that  diffusion  of  oxygen  into  and 
out  of  the  water  can  be  calculated.  Diffusion  was  negligible  in  all 
cases.  The  rate  of  oxygen  change  in  the  water  column  is  shown  in  the 
lower  part  of  the  figure.  The  area  under  the  segment  of  the  rate  curve 
which  is  bounded  by  the  solar  insolation  curve  (dashed  line)  represents 
daytime  net  photosynthesis.  The  area  under  the  other  rate  curve  repre- 
sents nighttime  respiration.  All  values  above  the  zero  line  indicate 


STATION  IB,  APRIL  30- MAY  I,  1977 


Figure  2.  Sample  set  of  curves  calculated 
from  measurements  of  dim'nal  changes  in 
oxygen  concentration 


a net  production  of  oxygen;  values  below  the  line  indicate  that 
utilization  of  oxygen  was  greater  than  production. 

11.  The  following  formula  was  used  to  estimate  gross  primary 
productivity: 
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gross  primary  productivity  = daytime  net  photosynthesis  + ^night 
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12.  Wetzel  indicates  that  the  relationship  between  day  and  night 
respiration  has  not  been  evaluated  satisfactorily  and  that  current 
methods  of  estimating  photosynthesis,  or  gross  primary  productivity,  at 
best  allow  only  general  comparisons  to  be  made. 

13.  The  amount  of  carbon  fixed  was  determined  by  assuming  a ratio 

of  2.67  to  1 for  oxygen  to  carbon  produced.  Light  and  dark  bottles 

provided  estimates  of  planktonic  productivity  and  respiration  in  the 

photic  zone  at  each  station.  These  values  were  subtracted  from  the 

values  obtained  in  the  diurnal  measurements  to  estimate  productivity  of 

the  macrophyte-epiphyte  complex  and  respiration  of  macrophytes, 

epiphytes,  and  benthos.  During  the  study  year,  Mitella  was  the  major 

species  of  macrophyte  in  the  lake.  Since  this  macroalga  does  not  have  a 

vascular  system,  the  use  of  diurnal  dissolved  oxygen  measurements  for 
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estimating  its  productivity  is  valid.  Hartman  and  Brown  have  pointed 
out  that  vascular  plants  store  much  of  the  oxygen  they  produce  within 
their  lacunar  system,  preventing  detection  by  the  diurnal  technique. 


Results 


lU.  The  annual  metabolic  pattern  shown  in  Figure  3 for  an  average 
square  metre  of  lake  surface  includes  three  peaks  of  gross  production: 
spring,  fall,  and  midwinter.  The  peaks, of  community  respiration  fol- 
lowed the  same  general  pattern,  but  rarely  did  the  rates  of  gross  produc- 
tion and  respiration  equal  each  other.  Net  community  production  (cal- 
culated by  subtracting  community  respiration  from  gross  production)  was 
jsitive  during  the  spring  and  late  summer  but  was  negative  at  all  other 
times.  Greatest  accumulation  of  organic  matter  occurred  in  the  late 
summer,  coinciding  with  the  large  pulse  of  gross  primary  productivity 
at  that  time.  A major  midsummer  decline  in  gross  primary  productivity 
as  well  as  in  respiration  also  occurred. 
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Figure  3.  Annual  pattern  of  production  and  respiration  calculated 
for  Lake  Conway.  (Measurements  were  made  monthly  at  nine  stations. 
Vertical  bars  represent  standard  errors  of  the  mean) 

15.  The  distribution  of  gross  productivity  for  different  depths 

and  times  is  shown  in  Figure  1*.  At  the  shallow  stations  (Figure  l*a), 

peaks  of  gross  productivity  occurred  in  the  spring  and  late  summer.  The 

larger  peak  was  in  the  spring  and  lasted  nearly  twice  as  long  as  the 

late  summer  peak.  The  maximum  spring  value  of  gross  productivity  at  the 

2 

shallow  depths  was  slightly  greater  than  2.0  g C/m  per  day. 

16.  At  the  medium-depth  stations  (Figure  Ub),  a spring  peak  of 

productivity  occurred,  but  there  was  no  late  summer  peak.  Highest 

productivity  at  these  stations  occurred  near  the  bottom,  and  maximum 
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values  were  slightly  greater  than  1.75  g C/m  per  day.  The  duration  of 
the  spring  productivity  pulse  was  about  the  same  as  that  for  the  shallow 
stations. 
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b.  Medium-depth  stations 
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c.  Deep  stations 

Figure  ^4.  Changes  in  distribution  of  productivity 
2 

(g  C/m  per  day)  with  depth 


17.  The  deep  stations  also  showed  a spring  pulse  of  productivity- 
near  the  bottom,  but  it  was  smaller  in  magnitude  and  shorter  in  duration 
than  the  others  (Figure  Uc).  Values  of  1.0  g C/m  per  day  were  observed 
several  times  throughout  the  year  at  depths  of  1 and  3 m. 

18.  Comparison  of  the  annual  metabolic  rates  of  the  three  commun- 
ity types  shows  that,  on  an  areal  basis,  the  pelagic  zone  had  the 
highest  gross  productivity  and  the  Vallisneria-Potamogeton  (shallow) 
zone  had  the  lowest  (Figure  5).  Communities  in  both  of  these  zones 


NCP  = Net  Community  Production 
CR  = Community  Respirotion 

NCP(»  34.3) 


VALLISNERIA-  HYDRILLA-  PELAGIC  WEIGHTED  MEAN 

POTAMOGETON  NIT  ELLA  ZONE  FOR  LAKE 

Figure  5-  Total  productivity  and  respiration  for 
shallow,  medium-depth,  and  deep  zones  in  Lake 
Conway  and  weighted  mean  for  entire  lake 

exhibited  positive  net  community  respiration  in  excess  of  gross  produc- 
tion. The  weighted  mean  metabolic  rates  for  the  entire  lake  show  a 
nearly  perfect  balance  between  gross  production  and  community  respira- 
tion.  The  gross  productivity  value  of  898  g C/m  per  year  is  in  the 
eutrophic  range  as  classified  by  Wetzel.1 


19.  If  the  areal  metabolism  values  are  divided  by  their  respec- 
tive depths,  an  estimate  of  the  metabolism  for  an  average  cubic  metre  in 
the  water  column  is  obtained  (Figure  6).  It  is  now  clear  that 


Figure  6.  Productivity  and  respiration  per  cubic  metre  for 
shallow,  medium-depth,  and  deep  zones  in  Lake  Conway 

productivity  on  a volumetric  basis  was  highest  in  the  shallow  zone  and 
lowest  in  the  pelagic  zone. 

20.  Temperature  measurements  during  1976  showed  that  deeper  areas 
of  the  lake  exhibited  a moderately  strong  stratification  which  began  in 
late  April  and  lasted  until  late  September  (Figure  7).  A much  shorter 
stratification  period  began  in  mid-February  and  lasted  until  late  March. 
Coinciding  with  the  longer  period  of  stratification  was  depletion  of 
dissolved  oxygen  in  hypolimnetic  waters  at  three  sites  which  were  deeper 
than  6.5  m.  This  anoxic  period  lasted  from  late  March  to  late  August. 

Discussion 


21.  Many  factors  may  be  interacting  to  produce  peaks  and  troughs 
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b.  Medium-depth  stations 


DEPTH,  (m) 


c.  Deep  stations 

Figure  7.  Changes  in  temperature  (°C)  over  time  and  depth 


of  productivity  and  respiration.  The  pattern  of  solar  radiation  on  Lake 
Conway  is  shown  in  Figure  8.  A notable  feature  in  this  pattern  is  the 


Figure  8.  Solar  radiation  of  Lake  Conway.  Solid  lines  represent 
monthly  averages , and  dashed  lines  represent  what  the  monthly 
averages  would  have  been  had  the  values  on  the  days  that  measure- 
ments were  made  been  typical.  (Data  provided  by  P.  Nawrocki, 

Florida  Solar  Energy  Center) 

midsummer  dip  in  sunlight  caused  by  the  daily  accumulation  of  thunder- 
heads  which  bring  late  afternoon  rains.  Estimated  inputs  of  total 
phosphorus  from  rainfall  and  runoff  are  shown  in  Figure  9.  Inputs  were 
highest  in  September,  when  20  percent  of  the  year's  rainfall  occurred, 
and  lowest  in  October  and  November,  when  only  3 percent  occurred. 

22.  The  midwinter  peaks  of  gross  productivity  and  respiration  do 
not  appear  to  be  related  to  any  one  variable.  However,  minimum  shading 
by  vascular  plants  and  lack  of  sufficient  sunlight  to  cause  photoinhibi- 
tion near  the  surface  may  have  brought  about  the  evenly  distributed, 
moderate  pulse  of  phytoplankton  productivity  seen  at  this  time.  Cold 
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Estimated  inputs  of  total  phosphorus  to  Lake  Conway 
(See  Table  6 for  assumptions) 


water  temperatures  in  February  probably  prevented  higher  levels  of 
productivity  and  respiration  from  developing. 

23-  The  trough  which  follows  this  midwinter  peak  may  have  been  a 
result  of  the  temperature  distribution  at  that  time.  As  the  lake  became 
slightly  stratified,  nutrients  in  the  photic  zone  may  have  been  depleted 
It  is  also  possible  that  photoinhibition  of  phytoplankton  may  have 
occurred  as  the  levels  of  solar  radiation  doubled  over  those  of  the 
previous  month. 

2b.  The  spring  peak  was  probably  due  to  increased  levels  of  light 
as  well  as  isothermal  conditions.  As  shown  previously,  gross  produc- 
tivity was  higher  near  the  bottom  for  all  depths  in  the  lake.  This 
suggests  that  the  major  portion  of  the  productivity  is  due  to  vascular 
plants  and  associated  epiphytic  algae,  especially  since  the  vascular 
plants  had  not  grown  to  the  surface  at  any  station. 

25-  The  midsummer  decline  in  the  lake's  metabolism  was  probably 
due  to  epilimnetic  nutrient  depletion  resulting  from  an  absence  of  free 
mixing  in  the  water  column.  The  slight  drop  in  sunlight  intensity 
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during  this  period  may  also  have  contributed  to  the  depressed  metabolism. 
It  should  be  noted  that,  since  the  hypolimnetic  waters  were  anoxic  during 
this  period,  the  diurnal  oxygen  technique,  which  measures  only  changes 
in  oxygen,  would  tend  to  underestimate  the  actual  respiration.  This 
problem  is  compensated  for  later  on,  however,  as  the  by-products  of 
anaerobic  metabolism  exert  a chemical  and  biological  oxygen  demand  as 
the  oxygen  is  replenished  to  the  hypolimnion  during  late  summer  turnover. 
Therefore,  part  of  the  community  respiration  values  seen  after  late 
summer  turnover  may  represent  a metabolic  time  lag. 

26.  The  late  summer  peak  of  productivity  can  also  be  related  to 
the  turnover  at  this  time  when  limiting  nutrients  such  as  phosphorus 
that  accumulate  in  the  anaerobic  hypolimnion  are  freely  mixed  within  the 
water  column.  The  depths  for  which  productivity  peaks  were  found  during 
this  time  suggest  that  the  phytoplankton  were  responsible  for  the 
majority  of  the  productivity. 

27.  A comparison  of  the  average  metabolic  rates  of  Lake  Conway 
with  those  of  other  lakes  of  varying  latitude  shows  no  clear  trend 
(Table  l).  Lake  morphology,  age,  and  basin  characteristics  are  of 
sufficient  importance  to  prevent  the  expression  of  a simple  latitudinal 
cline . 
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PART  III:  THE  LAKE  CONWAY  ECOSYSTEM  MODEL 
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28.  The  Lake  Conway  model  presented  in  this  report  is  based  on 
three  major  sources  of  information:  (a)  productivity  rates,  respiration 
rates,  and  relationships  among  submersed  plants  summarized  earlier  in 
this  report;  (b)  data  collected  by  other  researchers  participating  in 
the  Lake  Conway  baseline  studies;  and  (c)  estimates  and  assumptions 
obtained  from  an  extensive  literature  search.  Two  simulations  of  the 
ecosystem  are  presented.  The  first  is  a simulation  of  Lake  Conway  under 
"undisturbed"  conditions,  analyzing  conditions  that  were  measured  in  the 
lake  during  the  study  year  1976-1977  with  respect  to  the  variation  seen 
in  a 10-year  simulation.  The  second  is  a simulation  which  predicts  the 
effect  of  the  addition  of  white  amur  on  the  ecosystem  of  the  lake  over  a 
10-year  period.  The  model  itself  is  nonlinear  and  has  been  programmed 
in  CSMP  (Continuous  System  Modeling  Program),  a language  which  permits 
the  solution  of  difference  equations  on  a digital  computer.  The  program 
is  listed  in  Appendix  A. 

Data  and  Assumptions  Used  in  the  Model 

Plant  populations  in  Lake  Conway 

29.  Three  major  plant  components  are  included  in  the  model: 
phytoplankton,  epipelic  algae  (i.e.  that  living  upon  the  sediment),  and 
a component  which  includes  both  submersed  macrophytes  and  their  attached 
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epiphytic  algae.  Allen  has  shown  that  dissolved  organic  matter  re- 
leased by  mcarophytes  can  be  utilized  and  transformed  by  epiphytic  algal 
communities,  which  were  responsible  for  over  20  percent  of  the  total 
annual  production  in  a northern  lake.  The  macrophytes  and  their  as- 
sociated epiphytic  algae  therefore  appear  to  act  as  a complex  rather 
than  separately  and  have  been  considered  as  such  in  the  model. 

30.  All  macrophytic  plants  are  combined  into  one  population  unit 
in  the  model.  Since  productivity  data  were  not  obtainable  for  individ- 
ual species,  and,  at  present,  only  general  morphological  characteris- 
tics can  be  used  to  separate  the  species  by  function,  it  was  felt  that 
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separation  by  species  in  the  model  would  be  too  arbitrary  to  be  defen- 
sible. Moreover,  several  studies  on  the  white  amur  suggest  that  feeding 

preferences  of  the  fish  will  not  lead  to  the  exclusion  of  any  species. 

7 8 9 

Sills,  Michewicz,  Sutton,  and  Blackburn,  and  Fisch  and  Lyakhnovich, 

in  reviewing  studies  of  feeding  preferences,  conclude  that  the  diet  of 

the  white  amur  depends  almost  entirely  on  the  habitat  in  which  it  is 

found  and  that  it  will  eat  virtually  anything  in  the  absence  of  preferred 

food.  Prowse10  points  out  that  succulent  plants  with  low  fiber  content 

are  selected  preferentially.  During  the  study  year  at  Lake  Conway, 

Nitella  was  more  abundant  than  Hy dr ilia.  The  white  amur  is  known  to 
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readily  eat  Chara,  a closely  related  macroalga.  5 This  lack  of 
specificity  in  the  organism  of  interest  to  the  study  was  further  justi- 
fication for  not  distinguishing  between  plants  by  species  in  the  model. 

31.  Light  availability  for  photosynthesis  is  an  important  part  of 

the  equations  used  for  calculating  gross  primary  productivity  for  plant 

populations.  These  equations  focus  on  two  important  functions.  First, 

they  must  determine  the  amount  of  photosynt,.etically  active  radiation 

which  is  available  to  the  plant;  second,  they  must  predict  the  plant's 

photosynthetic  response  to  that  amount  of  light.  Several  authors  have 

derived  detailed  photosynthesis  equations  describing  still  or  mixed 
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phytoplankton  systems.  Few  photosynthesis  equations  for  systems 

containing  both  plankton  and  submersed  vascular  plants  have  been  pub- 
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lished.  One  such  equation,  in  the  model  proposed  by  Titus  et  al.,  is 
unsatisfactory  for  use  in  this  case  because  it  assumes  a constant  height 
for  the  submersed  vascular  plant  canopy.  In  Lake  Conway,  the  canopy 
height  changes  seasonally  and  so,  therefore,  does  the  amount  of  light 
available  to  it. 

32.  To  calculate  the  amount  of  light  available  to  the  vascular 
plants,  the  following  expression  is  used: 


Lq2  = S * (1  - R) 


e-(K2Z’  + K3Z'  • Q1/Z) 
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where 


S = solar  radiation  on  a horizontal  surface 

R = decimal  percentage  of  the  solar  radiation  reflected 

K ,K2  = constants  in  the  light  extinction  equation  for  Lake 
Conway  (which  is 


where 

I = radiation  at  depth  Z 

Lx 

Iq  = net  incident  radiation) 

Z'  = distance  from  the  water  surface  to  the  plant  canopy. 

This  varies  as  a function  of  cumulative  plant  biomass. 
The  function  was  determined  for  harvests  of  two  species 
of  plants  from  several  depths  at  five  stations.  The 
results  are  shown  in  Figure  10;  data  for  Potamogeton 
are  used  in  the  model 

K Z'  • Q^/Z  = empirical  relationship  which  reduces  the  light  avail- 
able to  the  submersed  plants  by  a function  of  the 
average  amount  of  phytoplankton  biomass  per  cubic 
metre  Qq/Z  and  the  depth  of  the  plant  canopy  Z'  . 

At  maximum  phytoplankton  levels , the  light  available 
for  the  vascular  plants  is  reduced  by  half 


33*  The  plant's  photosynthetic  response  to  the  available  light  is 
determined  by  the  plant's  efficiency  of  light  utilization  at  various 
light  levels.  The  predictive  efficiency  term  used  for  aquatic 
macrophytes  is 


JQ, 


1 + 


Lq  /L 
h2  ^2-50 


where 


En  = maximum  efficiency  for  naturally  occurring  light  conditions 

Lq  = predicted  amount  of  light  available  to  the  submersed 
2 plants 
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Figure  10.  Cumulative  biomass  of  two  species  of 
aquatic  plants  found  at  different  dept.'is 


L = light  level  at  which  the  plant  is  at  half  its  maximum 

^2-50  efficiency.  The  value  for  this  in  the  model  was  cal- 
culated from  photosynthesis  and  light  data  for  four  sub- 
mersed plant  speciesl? 

3^.  The  efficiency  term  is  multiplied  by  the  available  light  to 
give  production  in  kilocalories,  which  is  converted  to  grams  of  carbon 
by  a conversion  factor  of  9-3  kcal/g  C. 

35.  The  light  available  for  phytoplankton  is  predicted  by 


where 


constant 


self-shading  term  which  reduces  the  amount  of  available 
light  by  one  half  at  maximum  phytoplankton  biomass  levels 


,1 


s 


fgQg,  = function  describing  shading  of  phytoplankton  by  vascular 
plants.  At  maximum  vascular  plant  levels,  half  of  the 
light  available  to  phytoplankton  will  be  blocked 


3 6.  The  photo synthetic  response  of  phytoplankton  is  then  deter- 
mined by  solving  for  the  efficiency  of  light  use  with 


1 + lq  /lq 

^1  TL-50 


and  then  multiplying  by  available  light  to  give  productivity  in  kilo- 
calories and  converting  to  grams  of  carbon  by  a factor  of  9.3  kcal/g  C. 

l8 

E_  and  L_  were  calculated  using  data  from  Aruga. 

Ql  Ql-50 

37-  The  amount  of  light  available  for  epipelic  algae  is  predicted 
using  the  expression 


-(K  Z + K • Q ) 

Lq  = S • (1  - R)  • f3Q2  • e ' 

where 

f Q2  = function  describing  shading  of  epipelic  algae 
by  vascular  plants 

-(K2Z  + K • Qx) 

e = light  extinction  according  to  the  optical 

properties  of  the  water  as  well  as  shading  by 
phytoplankton 


38.  The  predictive  efficiency  term  for  epipelic  algae  is 


Q, 


1 + lq  /lq 

v j-50 


The  same  factor  for  conversion  to  grams  of  carbon  is  used. 

39.  The  average  relationship  between  net  photosynthesis  and 
temperature  in  three  species  of  aquatic  vascular  plants  is  shown  in 
Table  2.  The  original  values  were  converted  to  a relative  scale  in 
which  the  highest  value  for  the  three  reference  temperatures  became  one 

25 


and  the  other  two  values  were  fractions  of  that.  Data  summarized  in 
20 

Sculthorpe  suggest  that  Q (see  paragraph  63)  for  respiration  in 
aquatic  macrophytes  varies  from  1.32  to  3.1+8,  but  data  from  Ikusima 
suggest  that  for  hydrilla  may  be  at  the  lower  end  of  the  scale. 

Without  sufficient  data  for  all  species  of  plants  found  in  Lake  Conway, 
a value  of  1.4  was  selected.  The  values  calculated  from  this  relation- 
ship were  also  normalized,  and  the  normalized  values  are  the  ones  listed 
in  Table  2.  In  order  to  combine  temperature  relationships  for  net  pro- 
ductivity and  respiration  into  a factor  that  could  be  used  to  express 
temperature  dependence  of  gross  primary  productivity,  the  relative  values 
for  net  productivity  and  respiration  were  summed,  and  the  resulting 
value  for  gross  primary  productivity  was  in  turn  converted  to  a relative 
scale  (Table  2). 

1+0.  Temperature  relationships  for  phytoplankton  were  taken  from 

-j  ~|  ^ 

Wetzel's  interpretation  of  Aruga's  data.  The  same  relationships  were 
assumed  to  hold  for  epipelic  algae. 

1+1.  The  inputs  and  outputs  to  the  macrophyte  complex  are  outlined 
in  Figure  11.  Biomass  data  were  collected  by  researchers  from  the 
Florida  Department  of  Natural  Resources . Dry  weight  values  were  con- 

20 

verted  to  grams  of  carbon  by  assuming  0.1+0  g C per  gram  of  dry  weight. 
Total  biomass  of  submersed  aquatic  macrophytes  was  actually  about 
3l+  g C/m  in  Lake  Conway  in  January.  Odum  estimated  that  30  percent 
of  the  plant  biomass  in  Silver  Springs,  Fla.,  was  attributable  to 
epiphytic  algae;  in  the  model,  it  is  assumed  that  they  represent  25  per- 
cent of  the  biomass  of  the  complex.  The  biomass  of  the  macrophyte- 
epiphyte  component  is  therefore  estimated  to  be  1+5  g C/m  . Gross  pri- 
mary productivity  values  were  obtained  from  the  productivity  study. 

1+2.  Nutrients  were  not  considered  to  be  limiting  to  this  group  of 

organisms  since  most  of  the  macrophytes  have  root  systems.  Bristow  and 
23  2l+  25 

Whitcombe,  Schults  and  Malueg,  and  Denny  have  shown  that  uptake  of 
phosphorus  through  the  root  system  is  an  important  process  in  many 
rooted  submersed  plants.  Moreover,  Nichols  and  Keeney^  have  shown  that, 
although  Myriophyllum  spicatum  L.  can  take  up  nitrogen  through  its 
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Figure  11.  Inputs  and  outputs  affecting  biomass  of  aquatic 
macrophytes  and  associated  epiphytic  algae.  (Explanations 
of  terms  in  equations  are  listed  in  Appendix  B) 

leaves,  it  can  satisfy  all  its  requirements  through  root  uptake  from 
sediment . 

1*3.  Wetzel'1'  indicates  that  secretion  of  organic  carbon  by  macro- 
phytes may  vary  from  0.05  percent  to  over  100  percent  of  the  photo- 
synthate  but  that  most  values  lie  between  1 and  10  percent.  However, 
he  also  points  out  that  exchanges  of  inorganic  and  organic  matter  occur 
between  macrophytes  and  their  associated  epiphytic  algae.  In  the  model, 
it  is  assumed  that  no  organic  carbon  is  leached  into  the  surrounding 
water. 

1+i*.  Herbivory  is  estimated  to  constitute  only  5 percent  of  net 
primary  productivity  (8  percent  is  the  maximum  amount  estimated  by 


Wetzel'*')  and  much  of  this  is  presumed  to  be  concentrated  on  the  epiphytic 
algae.  Thus,  much  of  the  biomass  that  appears  as  net  primary  productiv- 
ity becomes  dead  organic  matter,  i.e.,  leaves  and  other  plant  parts  that 
die  in  situ.  In  the  model,  10  percent  of  the  biomass  each  month  is 
sloughed  off  as  dead  organic  matter,  and  an  additional  pulse  occurs  in 
early  September. 

U5.  Several  species  of  plants  in  Lake  Conway,  including  hydrilla, 
form  tubers.  These  are  included  as  a separate  state  variable  in  Fig- 
ure 11.  In  the  model,  these  tubers  are  assumed  to  start  forming  in 

October  and  to  begin  germination  in  late  April.  Haller,  Miller,  and 
27 

Garrard  have  found  that  tubers  of  hydrilla  form  in  the  bottom  sediment 

in  the  fall,  winter,  and  early  spring  months  and  germinate  in  the  late 

28 

spring,  summer,  and  early  fall.  Miller,  Garrard,  and  Haller  report 

2 

tuber  biomass  of  more  than  273  g/m  (fresh  weight)  in  a reservoir  in 

2 2 

Florida  or  approximately  11  g C/m  . A value  of  30  g C/m  per  month  was 
chosen  for  the  model  to  compensate  for  the  presence  of  tubers  of  other 
species  as  well. 

1+6.  The  phytoplankton  component  of  the  model  is  shown  in  Fig- 
ure 12.  Biomass  of  phytoplankton  in  January  was  calculated  to  be 
2 

0.8U  g c/m  for  Lake  Conway.  A team  from  the  Orange  County  Pollution 
Control  Department  provided  estimates  of  chlorophyll  a,  which  were  con- 
verted to  carbon  by  assuming  that  1 mg  of  chlorophyll  a is  equivalent  to 
50  mg  of  carbon.  Gross  primary  productivity  values  were  obtained  from 
the  productivity  study. 

1+7.  Zooplankton  and  fish  consume  phytoplankton  at  35  percent  of 

29 

net  primary  productivity.  Wetzel  et  al.  found  that  up  to  18  percent 
of  gross  primary  productivity  was  released  as  organic  carbon;  the 
annual  mean  percentage  was  7 percent.  Several  factors  affect  this 
release,  but,  in  their  study,  rates  seemed  to  be  highest  in  early  spring 
when  productivity  was  increasing  and  in  late  summer  when  it  began  de- 
creasing. A function  is  incorporated  in  the  model  to  allow  secretion  of 
18  percent  of  the  photosynthate  under  intermediate  gross  primary  produc- 
tivity levels.  This  decreases  gradually  to  1 percent  at  higher  and 

30 

lower  levels.  Golterman  estimates  that  20  percent  of  the  total  net 
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Figure  12.  Inputs  and  outputs  affecting  biomass  of  phytoplankton. 
(Explanations  of  terms  in  equations  are  listed  in  Appendix  B) 

productivity  may  reach  the  sediment;  sediment  deposition  accounts  for 
27  percent  in  the  model. 

48.  Data  on  epipelic  algae  were  not  available  for  Lake  Conway. 
However,  observations  during  the  study  indicate  that  this  is  an  impor- 
tant component  of  the  community.  The  relationships  that  were  derived 

31  32 

for  this  component  are  shown  in  Figure  13.  Hargrave  and  Gruendling 

2 

reported  gross  primary  productivity  levels  of  40  to  45  g C/m  per  year 
for  epipelic  algae  in  Marion  Lake,  British  Columbia.  Gross  primary 

2 

productivity  for  January  in  Lake  Conway  was  estimated  to  be  3.3  g C/m  . 

2 

The  January  standing  crop  was  estimated  to  be  1.0  g C/m  based  on  an 
arbitrary  turnover  time  of  10  days.  Nutrient  relationships  for  these 
algae  are  not  clear,  but  they  are  assumed  to  depend  on  dissolved 
orthophosphate  in  the  hypolimnion  because  of  the  depth  at  which  they 
grow.  Leaching  of  organic  carbon  was  assumed  to  follow  the  relationship 
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Figure  13.  Inputs  and  outputs  affecting  biomass  of  epipelic  algae. 
(Explanations  of  terms  in  equations  are  listed  in  Appendix  B) 

described  for  phytoplankton.  Consumption  is  Ul  percent  of  net  primary 
productivity. 

Animal  populations  in  Lake  Conway 

i*9.  Dynamics  of  zooplankton,  benthic  invertebrates,  and  four 
divisions  of  fish  (herbivorous,  young  and  adult  primary  predator,  and 
secondary  predator  fish)  are  included  in  the  model.  The  rates  of 
respiration,  consumption,  and  predation  of  these  organisms  are 

temperature-dependent.  These  relationships  are  incorporated  into  the 

33 

model  by  using  terms  reported  by  Patten,  Egloff,  and  Richardson  for 
organisms  in  an  Oklahoma  reservoir.  Although  the  effects  of  temperature 
on  consumption  and  digestion  are  treated  separately,  the  average  of  the 
two  values  is  used  in  the  model.  The  relationship  for  zooplankton  is 
Q equals  1.25;  the  following  terms  are  used  for  the  other  organisms: 


30 


Benthic  invertebrates 


Fish 


50.  The  values  calculated  with  these  terms  are  normalized  so  that 
full  activity  predicted  by  a term  in  an  equation  occurs  at  the  highest 
temperature.  The  fractions  of  full  activity  for  a year's  span  are  shown 
for  zooplankton,  benthic  invertebrates,  and  fish  in  Figure  lU.  The 
temperatures  used  were  recorded  during  the  diurnal  measurements  in  1976- 
1977  and  reflect  the  depths  most  commonly  occupied  by  the  organisms. 
Thus,  benthic  invertebrates  are  affected  only  by  temperatures  at  the 
bottom  of  the  lake,  while  the  zooplankton  and  fish  are  affected  by  the 
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Figure  lU.  Annual  variation  in  feeding  and  respiration  rates  of 
aquatic  organisms.  (Relationships  used  to  calculate  activity  are 
given  in  text.  Monthly  temperatures  used  in  these  relationships 
were  actual  measurements  from  different  depths  at  Lake  Conway) 


average  temperature  of  the  lake.  In  the  model,  the  rates  of  predation 
of  fish  and  benthic  invertebrates  are  programmed  as  functions  of 
temperature-dependent  percentages  of  activity  of  both  predators  and  prey. 
Zooplankton,  for  instance,  might  be  at  an  advantage  over  both  benthic 
invertebrates  and  fish  in  the  late  fall  and  early  spring  when  they  enjoy 
a higher  percentage  of  their  potential  activity.  Fish  are  dominant  over 
the  others  in  late  spring  through  midsummer  but  are  relatively  less 
active  by  midfall. 

51.  The  hypothesized  feeding  relationships  in  the  Lake  Conway 
ecosystem  are  shown  in  a matrix  in  Figure  15 . Initial  amounts  were 
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Figure  15.  Feeding  relationships  among  plants,  animals, 
and  detritus  in  the  Lake  Conway  ecosystem  model.  Numbers 
are  percentages.  The  percent  composition  of  an  item 
listed  along  the  top  in  the  diet  of  an  animal  listed  on 
the  left  side  is  given  in  the  lower  left  half  of  each 
square;  the  percent  of  each  item  listed  along  the  top  that 
is  consumed  by  the  animal  on  the  left  is  listed  in  the 
upper  right  half  of  each  square 
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1 


determined  arbitrarily,  and  many  of  these  were  readjusted  during  the 
initial  simulation  of  the  ecosystem.  The  type  of  plant  eaten  most 
universally  is  of  the  macrophyte-epiphyte  complex.  The  macrophytes  them- 
selves are  probably  only  lightly  grazed,  considerably  heavier  pressure 
being  directed  toward  the  epiphytic  algae.  This  group  is  grazed  by 
every  consumer  except  the  secondary  predator  fish  and  is  utilized  most 
heavily  by  zooplankton  and  benthic  invertebrates.  These  two  organisms 
are  the  most  universally  grazed  of  the  animals.  The  herbivorous  fish 
have  the  widest  food  spectrum  of  any  of  the  consumers,  with  a preference 
toward  epipelic  algae. 

52.  The  zooplankton  component  of  the  model  is  shown  in  Figure  16. 
Biomass  values  for  zooplankton  were  calculated  by  multiplying  numbers 
of  zooplankton  per  square  metre  by  the  dry  weight  of  each  individual  as 

34 

reported  in  Brezonik,  Morgan,  and  Shannon.  Total  dry  weight  for  each 
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Figure  l6.  Inputs  and  outputs  affecting  biomass  of  zooplankton. 
(Explanations  of  terms  in  equations  are  given  in  Appendix  B) 
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month's  sample  was  then  multiplied  hy  0.5  to  convert  the  value  to  grams 

2 

of  carbon.  Values  for  1976-1977  varied  from  0.15  to  0.W  g C/m  per 

2 

month.  This  is  lower  than  the  value  of  1.29  g C/m  reported  by  Nordlie 
35 

to  Grocki  for  Newman’s  Lake  in  Alachua  County,  Fla.  Annual  variation 
in  standing  crop  (expressed  as  a percentage  of  the  maximum  monthly  stand 
ing  crop)  is  shown  in  Figure  17- 


J FMAMJJASOND 
Time 

Figure  17.  Relative  variations  in  standing  crop  of 
zooplankton.  (Values  were  calculated  from  data  pro- 
vided by  a research  team  from  the  University  of 
Florida's  Department  of  Environmental  Engineering 
Sciences.  Possible  trend  of  values  within  data  gap 
are  represented  by  the  dashed  line) 

53.  Several  studies  summarized  by  Wetzel1  show  that,  in  general, 
10  to  75  percen  of  the  energy  ingested  by  zooplankton  is  assimilated. 
Gulati"5  estimated  36  percent  efficiency  for  herbivorous  zooplankton  in 
Dutch  lakes.  Assimilation  in  the  Lake  Conway  model  is  estimated  to  be 

3k 


47  percent  of  ingestion.  Gulati  also  estimated  that  respiration  is 

20  to  30  percent  of  the  amount  of  energy  ingested,  although  Wetzel 

37 

listed  values  varying  from  27  to  93  percent.  McAllister  estimated 
that  respiration  is  roughly  10  percent  of  the  standing  crop  per  day.  If 
respiration  is  assumed  to  account  for  two  thirds  of  the  energy  assim- 
ilated (or  31  percent  of  energy  ingested),  the  respiration  rate  in  the 
model  is  12  percent  of  the  standing  crop  per  day.  The  general  distri- 
bution of  energy  resources  is  therefore  in  agreement  with  those  from 
these  previous  studies.  Wetzel^"  summarized  several  zooplankton  studies 
which  showed  turnover  (net  annual  productivity  divided  by  biomass)  vary- 
ing between  9 and  25  per  year  in  eutrophic  Polish  lakes,  with  values  of 
l4  and  15  per  year  for  other  studies.  Turnover  in  the  Lake  Conway  popu- 
lation, using  the  assumptions  given  and  the  January  biomass  estimate,  is 

21  per  year,  which  is  reasonable  when  climatic  differences  are  taken 
into  account. 

54.  The  benthic  invertebrate  component  of  the  model  is  shown  in 
Figure  18.  Four  months  of  data  (April,  May,  July,  and  September)  were 
available  for  calculation  of  benthic  invertebrate  biomass.  Biomass  estl 

i 38  39 

mates  summarized  in  Wetzel  from  Junk'3  and  Sapkarev  and  Tocko0  were 

used  to  convert  numbers  of  organisms  to  dry  weights  of  sample.  Dry 
weights  were  then  converted  to  grams  of  carbon  by  multiplying  by  0.5. 
Biomass  varied  from  1.35  to  2.27  g C/m  . Grocki  estimated  1.45  g C/m 
for  Lake  Kissimmee,  while  citing  data  from  other  authors  which  varied 
between  1.07  and  3.00  for  different  parts  of  the  lake.  Since  no  sea- 
sonal trends  were  apparent  from  the  Lake  Conway  data  (the  lowest  value 

was  for  September  and  the  highest  for  April),  an  intermediate  value  of 
2 

1-73  g C/m  was  selected  to  represent  initial  conditions  in  January. 

40 

55.  In  a review  of  feeding  studies  of  aquatic  insects,  Cummins 

35 

stresses  the  heavy  dependence  of  this  group  on  detritus.  Grocki  esti- 
mated, after  reviewing  several  reports,  that  algae  comprised  30  percent 
and  detritus  JO  percent  of  the  diet  of  benthic  invertebrates  in  Lake 

Kissimmee.  It  is  assumed  in  the  Lake  Conway  model  that  all  the  inverte- 

4l 

brates  are  detrivorous  except  Chaoborus , which  preys  on  zooplankton, 

42 

and  snails,  which  graze  the  epiphytes  on  macrophytic  plants. 
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Figure  l8.  Inputs  and  outputs  affecting  biomass  of  benthic 
invertebrates.  (Explanations  of  terms  in  equations  are 
given  in  Appendix  B) 


Proportions  of  the  overall  diet  comparable  to  the  biomass  of  these 
organisms  are  specified  in  the  model.  Detritus  comprises  68  percent  of 
the  diet. 

56.  Ingestion  rates  summarized  by  Cummins^  indicate  that  2 to 

23  percent  of  the  body  weight  of  aquatic  insects  is  ingested  per  day. 

The  rate  used  in  the  model  is  2 percent.  Assimilation  efficiency  is 
. 1*3 

69  percent,  a value  determined  by  Zimmerman,  Wissing,  and  Rutter  for 
mayfly  larvae. 

57-  Fish  were  collected  by  researchers  from  the  Florida  Game  and 
Freshwater  Fish  Commission.  Data  from  block  net  samples  taken  in  Lake 
Conway  in  June  were  used  in  the  model.  The  fish  found  in  Lake  Conway 
are  listed  in  Tables  3-5.  Assignment  to  trophic  group  (herbivore,  pri- 
mary predator,  or  secondary  predator)  was  based  on  the  information  given 


in  these  tables.  Since  young  primary  predator  fish  eat  a substantially 
different  diet  from  the  adults,  they  are  separated  as  another  state 
variable  and  assumed  to  eat  the  same  diet  items  as  the  herbivorous  fish. 
Relative  proportions  of  the  diet  items  in  each  group  were  assigned 
before  simulations,  but  many  were  changed  during  simulation. 

58.  Total  fish  biomass  in  Lake  Conway  is  estimated  to  be 

2 65 

1-39  g C/m  by  assuming  1 g C per  10  g of  dry  weight.  This  biomass 

33 

value  is  approximately  half  of  what  Patten,  Egloff,  and  Richardson 

33 

estimated  for  an  Oklahoma  reservoir.  However,  Grocki  estimated  only 

2 2 
l6  g C/m  for  the  littoral  zone  of  Lake  Kissimmee  and  2 g C/m  for  the 

66 

limnetic  zone  using  data  from  Wegener,  Williams,  and  Holcomb.  Hcrrel* 

2 

estimated  0.7  to  3.5  g C/m  for  Lake  Trafford,  Fla.  Assimilation  effi- 
ciencies for  these  groups  of  fish  were  determined  by  the  following  proce- 
68 

dure.  Winberg  estimates  that  80  percent  of  the  food  consumed  by  wild 

69 

fish  is  metabolizable,  and  Balogh  has  pointed  out  that  nonselective 
feeders  have  lower  assimilation  efficiencies  than  more  highly  selective 
consumers.  In  the  Lake  Conway  model,  assimilation  efficiency  has  been 
set  at  O.69  for  the  least  selective  group,  the  herbivorous  fish,  at  0.80 
for  the  primary  predator  fish,  and  at  1.0  for  the  secondary  predator 
fish. 

59*  The  dynamics  of  the  herbivorous  fish  component  are  shown  in 

2 

Figure  19.  Biomass  is  estimated  to  be  0.07  g C/m  or  5 percent  of  the 

total  biomass  of  fish.  Using  an  arbitrary  turnover  time  of  4 years, 

2 

net  productivity  of  this  group  of  fish  is  estimated  at  0.0015  g C/m  per 

month.  For  all  fish,  removal  by  fishing  is  ignored,  and  predation  is 

considered  to  be  the  only  significant  cause  of  death  in  herbivorous  and 

70 

primary  predator  fish.  In  the  model,  the  herbivorous  fish  depends  most 
heavily  on  epipelic  algae  but  equally  on  all  other  food  items.  The  feed- 
ing rate  is  1 percent  of  body  weight  per  day  in  January;  Patten,  Egloff, 
33 

and  Richardson  estimated  no  more  than  10  percent  per  day. 

60.  Primary  predator  fish  are  divided  into  two  groups,  as  shown 
in  Figure  20.  The  young  fish  are  born  in  April,  and  most  have  advanced 


# 
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Unpublished  report  cited  by  Carlson  and  Duever. 


into  adult  age  classes  by  the  end  of  the  year.  As  young,  they  eat  the 
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same  food  items  as  the  herbivorous  fish.  LeCren  observed  that  perch 
stop  feeding  on  plankton  after  they  have  grown  to  100  mm  in  length,  so 
biomass  of  young  fish  in  Lake  Conway  was  determined  by  adding  the 
weights  of  all  fish  less  than  100  mm  except  the  minnows,  killifish,  and 
silverside.  Seventeen  percent  of  the  primary  predator  fish  fell  into 
this  category. 

61.  Primary  predator  fish  comprise  6k  percent  of  the  total  fish 

population.  Turnover  time  for  bluegills  in  northern  Indiana  is  slightly 
72 

more  than  1 year.  Turnover  time  of  the  adult  primary  predator  fish 
in  the  model  is  estimated  to  be  about  2 years.  Primary  predator  fish 
feed  most  heavily  on  the  macrophyte-epiphyte  complex,  moderately  on 
benthic  invertebrates,  and  less  on  epipelic  algae  and  zooplankton. 

62.  Secondary  predator  fish  (Figure  21 ) comprise  31  percent  of 
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Figure  21.  Inputs  and  outputs  affecting  biomass  of 
secondary  predator  fish.  (Explanations  of  terms  in 
equations  are  given  in  Appendix  B) 


the  total  fish  biomass.  Turnover  time  is  estimated  to  be  U years. 
These  fish  depend  most  heavily  on  primary  predator  fish  and  moderately 
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on  benthic  invertebrates;  only  2 percent  of  their  diet  is  derived  from 
the  herbivorous  fish. 

Detritus 


I 

J 


6 3.  All  dead  material  sinks  into  the  detritus,  from  which  respira- 
tion and  consumption  by  benthic  invertebrates,  herbivorous  fish,  and 
young  primary  predator  fish  are  the  only  outflows.  The  dynamics  of 
detritus  flows  are  shown  in  Figure  22.  The  variable  rates  of  loss 
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Figure  22.  Inputs  and  outputs  affecting  amount  of  detritus. 
(Explanations  of  terms  in  equations  are  given  in  Appendix  B) 


of  respiratory  products  are  discussed  in  the  following  section. 
Phosphorus  dynamics  in  Lake  Conway 

6b.  Levels  of  dissolved  orthophosphate  in  the  epilimnion  and  the 
hypolimnion  are  included  in  the  model  as  state  variables.  Data  on  phos- 
phorus levels  in  Lake  Conway  were  monitored  by  researchers  from  the 


1*0 


Orange  County  Pollution  Control  Department.  The  phosphorus  inputs  to 
Lake  Conway  through  rainfall  and  runoff  are  shown  in  Table  6.  Ortho- 
phosphate is  estimated  to  be  39  percent  of  total  phosphorus,  based  on 
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actual  measurements.*  Hendry  and  Brezonik  reported  3 6 percent  for  an 
Alachua  County  site  during  1976.  All  flows  of  phosphorus  in  this  model 
are  of  phosphorus  as  orthophosphate. 

65.  Dynamics  of  phosphorus  flow  in  the  epilimnion  are  shown  in 
Figure  23.  This  component  is  affected  by  the  input  of  phosphorus  from 
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Figure  23.  Inputs  and  outputs  affecting  level  of  orthophos 
phate  in  epilimnion.  (Explanations  of  terms  in  equations 
are  given  in  Appendix  B) 


rainfall  and  runoff,  by  mixture  with  the  hypolimnion  in  the  late  summer 
or  fall  and  spring  during  turnover,  and  by  uptake  and  excretion  of  phos 
phorus  by  plants  and  animals.  However,  few  data  could  be  found  on  the 
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magnitudes  and  controls  of  these  amounts.  Consequently,  the  flow  of 
orthophosphate  was  at  first  assumed  to  be  proportional  to  the  flow  of 
carbon,  with  the  C:P  ratio  of  roughly  100:1  used  for  conversion.  In 
many  cases,  this  had  to  be  altered  in  order  to  derive  a balanced  phos- 
phorus budget.  These  cases  are  noted  below. 

66.  Release  of  phosphorus  may  occur  in  three  ways:  passive 
release  as  a normal  accompaniment  to  plant  respiration  and  to  egestion 
and  excretion  in  animals,  active  secretion  by  plants  along  with  organic 
carbon  and  associated  compounds;  and  leaching  from  dead  materials. 

67.  The  pumping  of  phosphorus  from  the  sediment  to  the  water  by 
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macrophytes  has  been  documented  in  several  instances.  As  mentioned 

earlier,  the  secretion  of  organic  compounds  is  not  believed  to  be  de- 
tectable when  a community  of  epiphytic  algae  is  present  on  the  plants. 
Losses  with  respiration  are  considered  to  release  only  half  the  propor- 
tional levels  of  phosphorus  in  plants.  Mineralization  from  dead  leaves 

is  probably  a relatively  slow  process,  although  there  are  few  data  to 
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document  this.  Pieczynska  indicates  that  two  species  of  Potamogeton 
required  7 to  lU  days  to  lose  6 to  92  or  95  percent  of  their  biomass. 

In  the  model,  it  is  assumed  that  only  5 percent  of  the  phosphorus  in 
dead  material  leaches  into  the  epilimnion  before  the  material  sinks  to 
the  bottom. 

68.  Phosphorus  is  released  from  phytoplankton  through  both  res- 
piration and  secretion  of  organic  matter.  Only  half  the  normal  comple- 
ment of  phosphorus  is  judged  to  be  released  through  respiration.  Phos- 
phorus uptake  is  considered  to  be  proportional  to  gross  primary 
productivity  and  phosphorus  secretion  to  the  active  release  of  organic 
compounds. 

69.  For  zooplankton,  25  percent  of  the  proportional  amount  of 

phosphorus  is  accounted  for  in  egestion  and  excretion.  Kitchell,  Koonce, 
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and  Tennis  report  that  bluegill  remineralize  less  than  1 percent.  The 
value  was  therefore  set  at  0.5  percent  of  egested  and  excreted  materials 
for  all  fish.  The  contribution  of  zooplankton  is  substantially  more 
than  this. 

70.  It  is  assumed  that  complete  turnover  between  epilimnion  and 


hypolimnion  will  be  achieved  in  1 month.  The  two  compartments  are  begun 
at  the  same  value,  since  measurements  by  the  Orange  County  Pollution  Con 
trol  Department  showed  no  consistent  differences. 

71.  Phosphorus  flows  in  the  hypolimnion  are  shown  in  Figure  2k. 
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Figure  2k.  Inputs  and  outputs  affecting  level  of  orthophosphate 
in  hypolimnion.  (Explanations  of  terms  in  equations  are  given 

in  Appendix  B) 


Release  of  phosphorus  from  the  sediments  involves  two  pathways.  It  is 

believed  to  be  fastest  when  the  lake  is  anaerobic  during  the  summer. 

8l 

This  assumption  is  based  on  studies  in  the  Great  Lakes  where  most 
of  the  phosphorus  is  commonly  adsorbed  onto  clays  and  ferric 
hydroxides.  If  the  lake  is  well-oxygenated,  the  top  few  millimetres 
of  the  sediments  are  oxidized,  and  iron,  manganese,  and  phosphate  are 
trapped.  The  phosphorus  adsorbs  onto  and  complexes  with  ferric 
oxides  and  hydroxide.  If  this  zone  remains  oxidized,  nothing 


escapes  from  the  sediments.  However,  when  oxygen  levels  at  the  inter- 
face between  water  and  sediments  fall  below  1 mg/X,,  orthophosphate  ions 
are  transferred  from  the  sediment  to  the  water.  Release  during  aerobic 

periods  will  be  much  lower,  but  activity  by  benthic  invertebrates  will 
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allow  some  phosphorus  to  be  released.  Stumm  and  Leckie  estimated  the 

2 

release  rate  for  undisturbed  sediment  in  the  summer  to  be  1 g/m  per 
year.  The  rate  in  the  model  is  assumed  to  be  150  percent  of  this, 
because  of  higher  temperatures,  and  the  release  rate  in  aerobic  waters 
is  arbitrarily  assumed  to  be  one  third  of  the  anaerobic  rate. 

72.  Half  of  the  phosphorus  leaching  from  dead  phytoplankton  sink- 
ing from  the  epilimnion  is  released  into  the  hypolimnion;  5 percent  of 
the  phosphorus  in  dead  organic  material  sloughed  from  macrophytes  is 
also  released.  The  phosphorus  flows  involving  epipelic  algae  in  the 
hypolimnion  derive  from  the  same  assumptions  that  governed  the  phyto- 
plankton in  the  epilimnion.  Like  the  phytoplankton,  half  the  phosphorus 
in  the  dead  material  leaches  into  the  hypolimnion  before  becoming 
detritus . 

73.  The  phosphorus  from  egested  and  excreted  products  by  the  ben- 
thic invertebrates  is  released  into  the  hypolimnion.  Fish  bones  and 
scales  mineralize  very  slowly;  50  percent  of  the  phosphorus  in  fish 
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biomass  is  in  this  form.  ’ Secondary  predator  fish  are  the  only 
group  making  a substantial  contribution  to  this  pathway. 


Simulation  of  the  Ecosystem 

jk.  The  Lake  Conway  ecosystem  was  simulated  for  a 10-year  period 
using  January  measurements  for  initial  conditions.  Although  it  is  un- 
likely that  Lake  Conway  is  currently  in  steady  state,  the  goal  of  the 
modeling  effort  is  the  achievement  of  a steady  state  simulation  for  a 
10-year  period.  Without  knowing  the  histories  of  each  component,  it 
must  be  assumed  that  the  correlations  observed  among  components  during 
the  study  year  can  represent  a stable  condition. 

75*  Simulated  changes  in  biomass  of  the  macrophyte-epiphyte 
complex  are  shown  in  Figure  25.  Summer  biomass  in  the  tenth  year  is 
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Figure  25.  Simulated  pattern  of  change  in 
macrophyte-epiphyte  complex  biomass  under 
undisturbed  conditions  for  10  years 

12  percent  greater  than  in  the  first  year,  but  the  annual  pattern  of 
growth  and  decline  is  stable.  The  simulated  changes  in  the  second  year 
are  shown  in  Figure  26,  along  with  seven  measurements  reported  by 
researchers  from  the  Florida  Department  of  Natural  Resources.  The  simu- 
lation does  not  show  biomass  in  the  summer  to  reach  nearly  the  extreme 
suggested  by  the  September  data  point.  Nor  does  the  simulation  show  the 
secondary  biomass  increase  in  December.  The  significance  of  the  higher 
September  biomass  will  be  evaluated  when  the  summer  data  are  available; 
the  small  peak  in  December  is  probably  not  significantly  different  from 
the  other  fall  and  winter  points.  Further  information  on  rates  of 
biomass  sloughing  by  the  macrophytes,  nutrient  uptake  mechanisms,  and 
photosynthetic  responses  to  variable  light  conditions  might  provide 
clues  to  discrepancies  between  assumptions  and  reality. 

j6.  The  simulated  changes  in  phytoplankton  biomass  are  shown 
in  Figure  27-  One  year's  data  from  the  Orange  County  Pollution  Con- 
trol Department  and  a partial  set  of  data  from  researchers  from  the 
University  of  Florida's  Department  of  Environmental  Engineering 
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Figure  2J.  Simulated  changes  in  phytoplankton  bio- 
mass under  undisturbed  conditions  for  10  years 


Month 

Figure  28.  Simulated  changes  in  phytoplankton  biomass  for 
1 year.  (Biomass  estimates  made  from  data  reported  by  teams 
from  the  Orange  County  Pollution  Control  Department  and  from 
the  University  of  Florida's  Department  of  Environmental 
Engineering  Sciences  are  also  shown) 

Sciences*  are  shown  in  Figure  28.  The  latter  set  of  da^a  was  collected 
under  a more  restricted  sampling  scheme.  Plotted  with  these  curves 
are  biomass  values  from  the  fourth  year  of  the  simulation.  The  simula- 
tion underestimates  summer  biomass  predicted  by  both  sets  of  data, 
although  it  follows  the  pattern  of  increase  shown  in  the  data  from  the 
Orange  County  Pollution  Control  Department. 

77.  Nutrient  relationships  between  the  phytoplankton  and  the 
surrounding  water  are  a sensitive  area  in  the  model,  and  the  factors 
controlling  rates  of  uptake  and  loss  are  not  clearly  defined  in  the 
literature.  In  further  simulations,  the  effect  of  nutrients  will  be 
modeled  as  a limiting  factor;  in  the  model  at  present,  only  availability 

* J.  Fox  et  al. , personal  communications. 
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of  phosphorus  controls  phytoplankton  growth.  There  Is  little  basis,' 
however,  for  accurately  determining  the  dynamics  of  this  relationship. 
Seasonal  changes  that  have  not  been  accounted  for  in  the  relationship 
between  chlorophyll  a and  biomass  may  also  affect  the  biomass  estimate 
itself. 

78.  Winter  biomass  of  epipelic  algae  evolves  to  be  about  50  per- 
cent higher  than  predicted  in  the  initial  conditions  (Figure  29). 


Figure  29-  Simulated  changes  in  biomass  of 
epipelic  algae  under  undisturbed  conditions 
for  10  years 

Interrelationships  between  this  component  and  nutrients  in  the  bottom 
sediment,  water,  and  macrophytic  plants  (epipelic  algae  frequently 
grow  thickly  around  the  plants)  need  better  definition.  However,  little 
work  has  been  published  on  epipelic  algae,  and  measurements  of  its 
biomass  in  Lake  Conway  are  not  available. 

79.  Neither  zooplankton  nor  benthic  invertebrates  display  a 
stable  simulation  pattern  until  the  fifth  year  of  simulation.  Zooplank- 
ton (Figure  30)  do  not  show  the  midsummer  decrease  in  biomass  recorded 
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Figure  30.  Simulated  changes  in  biomass  of  the 
zooplankton  population  under  undisturbed  condi- 
tions for  10  years 


in  Lake  Conway  (Figure  IT).  A biomodal  pattern  such  as  apparently 
occurs  in  Lake  Conway  is  common  in  northern  temperate  lakes,  where  many 


of  the  smaller  species  are  subject  to  predation  by  larger  zooplankton 
in  the  summer  while  other  species  aestivate.  The  importance  of 


fish  predation  may  be  underestimated  in  the  model;  several  studies  cited 
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by  Wetzel  stress  the  importance  of  fish  predation  on  size  distribution 
among  zooplankton. 

80.  The  small  seasonal  changes  in  biomass  shown  by  the  simulation 

U0 


of  benthic  invertebrates  (Figure  31)  are  probably  reasonable.  Cummins 
points  out  that  temperature  is  the  primary  control  over  feeding  and 
respiration  rates  when  food  sources  are  both  abundant  and  consistent  all 
year.  Four  measurements  of  biomass  of  benthic  invertebrates  in  Lake 


Conway  (taken  by  researchers  from  the  University  of  Florida's  Department 


of  Environmental  Engineering  Sciences)  varied  from  1.35  g C/m  in 


September  to  2.27  g C/m  in  April,  with  intermediate  values  of  1.53  in 
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BENTHIC  INVERTEBRATES 


Figure  31.  Simulated  changes  in  biomass  of  the 
benthic  invertebrate  population  under  undis- 
turbed conditions  for  10  years 

May  and  1.76  in  July.  One  sample  was  taken  in  each  of  these  months,  and 
the  number  of  individuals  in  the  samples  ranged  from  8 to  23b.  The 
variation  is  therefore  probably  not  significant.  This  simulated  pattern 
is  not  satisfactory,  however,  and  seems  to  be  a product  of  the  insta- 
bility that  occurs  early  in  the  simulation  of  the  zooplankton  component. 
Assumptions  of  productivity  and  respiration  rates  for  this  group  need  to 
be  reexamined  also  since  the  turnover  rate  does  not  fail  within  the 
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range  determined  by  Waters  in  a review  of  several  studies. 

8l.  Data  on  fish  populations  were  available  for  only  1 month,  so 
it  is  difficult  to  assess  the  results  of  the  simulations.  Changes  in 
biomass  of  all  fish  are  shown  in  Figure  32.  Simulations  of  herbivorous, 
primary  predator,  and  secondary  predator  fish  are  shown  in  Figures  33-35, 
respectively.  Levels  of  both  primary  predator  fish  and  herbivorous  fish 
vary  seasonally  and  rise  above  the  initial  levels  used  in  the  model. 
Irregularities  in  the  simulations  of  these  fish  appear  to  be  related  to 


Figure  32.  Simulated  changes  in  biomass  of  the 
fish  population  under  undisturbed  conditions 
for  10  years 


HERBIVOROUS  FISH 


Figure  33.  Simulated  changes  in  biomass  of  the 
herbivorous  fish  population  under  undisturbed 
conditions  for  10  years 


PRIMARY  PREDATOR  PISH 


Figure  34.  Simulated  changes  in  biomass  of  the 
primary  predator  fish  population  (both  young 
and  adult)  Tinder  undisturbed  conditions  for 
10  years 


Figure  35.  Simulated  changes  in  biomass 
of  the  secondary  predator  fish  popula- 
tion under  undisturbed  conditions  for 
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the  early  instability  of  the  zooplankton  and  benthic  invertebrate  popula- 
tions. Seasonal  variations  are  damped  in  the  secondary  predator  fish, 
and  this  group  does  not  show  an  increase  in  biomass. 

82.  Simulations  of  changes  in  orthophosphate  in  the  epilimnion 
and  in  the  hypolimnion  are  shown  in  Figures  36  and  37,  respectively. 

No  data  are  available  for  verification  of  these  simulations,  since 
concentrations  during  most  of  the  simulation  remain  below  the  limits  of 
detection  of  the  analyzer  used  by  the  Orange  County  Pollution  Control 
Department.  Levels  never  rose  above  this  concentration  in  Lake  Conway 
during  the  study  year.  The  simulations  show  decreases  in  both  the 
epilimnion  and  the  hypolimnion  in  the  spring  and  a second  decrease  in 
the  hypolimnion  in  late  summer.  Despite  the  fact  that  mixing  is  pro- 
grammed to  occur  from  late  summer  to  midwinter  (one  complete  turnover 
occurs  each  month  during  this  period),  the  two  layers  maintain 
distinct  patterns.  Detritus  (Figure  38)  doubled  over  the  10-year  period. 

83.  The  present  model  is  an  approximation  of  ecological  relation- 
ships in  Lake  Conway,  but  does  not  yet  represent  these  conditions  as 
closely  as  is  desired.  Zooplankton  are  the  only  group  whose  simulated 
pattern  of  seasonal  changes  in  biomass  is  unrealistic;  complete  data 

for  several  other  groups  were  not  available  when  this  report  was  written, 
however.  Although  simulated  patterns  of  change  in  macrophytes  and  in 
phytoplankton  are  consistent,  the  annual  productivity  and  respiration 
rates  rise  considerably  higher  than  the  levels  measured  and  programmed 
for  the  initial  year.  This  discrepancy  affects  both  uptake  and  release 
of  phosphorus  and  may  account  for  the  increases  seen  in  epipelic  algae 
and  in  all  the  animals  except  the  secondary  predator  fish. 

84.  The  importance  of  the  epipelic  algae  to  the  stability  of  the 
model  appeared  during  earlier  simulations,  when  phosphorus  released  by 
the  epipelic  algae  during  respiration  was  inadvertently  doubled.  This 
led  to  a five- fold  increase  in  herbivorous  fish  biomass,  a decrease  in 
benthic  invertebrate  biomass  (to  half  the  level  reached  in  the  simula- 
tion discussed  above),  and  a dramatic  decrease  in  secondary  predator 
fish.  When  the  model  has  been  stabilized  to  produce  a 10-year  steady 
state  simulation,  a structured  set  of  sensitivity  analyses  will  be 
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ORTHOPHOSPHATE  IN  EPtUMNION 
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Figure  36.  Simulated  changes  in  levels  of  ortho- 
phosphate in  the  epilimnion  under  undisturbed  con 
ditions  for  10  years 


Figure  37-  Simulated  changes  in  levels  of  ortho- 
phosphate in  the  hypolimnion  under  undisturbed 
conditions  for  10  years 
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Figure  38.  Simulated  changes  in  amount  of 
detritus  in  Lake  Conway  under  undisturbed 
conditions  for  10  years 


conducted  to  evaluate  the  effect  of  changes  in  each  component  on  the 
stability  of  other  components  and  of  the  entire  model. 

t I 
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PART  IV:  SIMULATED  INTRODUTION  OF  WHITE  AMUR 

Dynamics  of  White  Amur  Population 

85.  Much  of  the  published  information  on  the  white  amur  exists  in 

Russian  literature.  Many  of  the  studies  that  have  been  done  on  basic 

biological  features  of  the  white  amur  have  been  translated  and  summarized 

9 

by  Fischer  and  Lyakhnovich,  whose  interpretation  is  used  in  this 
analysis. 

86.  Growth  efficiency  data  for  white  amur  are  shown  in  Table  7- 
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The  first  set  of  data  summarizes  the  opinion  of  Fischer  and  Lyakhnovich 

that  three  distinct  growth  phases  are  discernible  in  white  amur:  during 

the  first  2 years,  fish  can  double  their  weight  within  a year;  from  the 

second  to  fifth  year,  growth  efficiency  is  3 percent,  and  body  weight 

may  increase  as  much  as  30  percent  per  year;  and  after  the  fifth  or 

sixth  year,  when  sexual  maturity  has  been  reached,  only  about  1 percent 

of  the  food  intake  is  used  for  body  growth.  Variability  in  the  second 

set  of  data  listed  is  due  primarily  to  different  weeds  which  were  fed 

88 

the  fish  to  determine  weed-specific  growth  efficiencies. 

87.  Weight-specific  digestion  rates  are  listed  in  Table  8. 

9 

Fischer  and  Lyakhnovich  summarize  data  showing  that  daily  growth  in- 
creases vary  from  29  percent  of  initial  body  weight  in  fry  to  2 percent 
in  young  fish  and  0.07  percent  in  mature  fish. 

88.  The  white  amur  component  of  the  model  is  shown  in  Figure  39- 

The  values  are  percentages  of  ingested  food,  assuming  that  the  fish  eat 

. 9 

100  percent  of  their  body  weight  daily.  Fischer  and  Lyakhnovich  report 
assimilation  efficiencies  ranging  from  31.2  percent  of  ingestion  to 
90.1  percent;  an  average  of  60.65  percent  is  used  in  the  model.  They 
also  estimate  that  85.5  percent  of  assimilated  food  is  lost  as  respira- 
tion or  death.  The  age-varying  growth  rates  are  incorporated  into  the 
model  by  adjusting  respiratory  losses  from  85.5  percent  for  a young  fish 

to  100  percent  for  a 5-  to  6-year-old  fish,  or  one  which  weighs  10  kg. 

9 

Fischer  and  Lyakhnovich  indicate  that  fish  can  reach  a size  of  32  kg, 
but  that  leveling  off  generally  occurs  closer  to  10  kg. 


56 


■ K ■ K.4,‘  W T.  ' K, 44 ' W K«  ■ ^ ' 0B 


• - K 


Figure  39-  Inputs  and  outputs  affecting  biomass  and 
number  of  white  amur.  (Egestion  and  assimilation  are 
expressed  as  percentages  of  food  ingested,  respira- 
tion is  expressed  as  percentage  of  food  assimilated. 

Mortality  is  expressed  as  the  daily  rate  of  loss  of 
numbers  of  fish) 

89.  Feces  and  dead  white  amur  are  assumed  to  sink  to  the  sediment 
with  some  remineralization  of  the  latter  releasing  into  the  hypolimnion. 


Simulated  Effects  of  White  Amur 


90.  The  introduction  of  7000  white  amur,  estimated  to  weigh 
1+50  g each,  was  simulated  with  the  Lake  Conway  ecosystem  model.  Changes 
in  the  number  and  biomass  of  the  fish  are  shown  in  Figures  1*0  and  1*1, 
respectively.  The  death  rate  of  the  fish  is  constant,  reducing  them  to 
about  1+3  percent  of  their  initial  level  in  10  years.  Biomass  increases 
rapidly,  however,  until  the  second  year.  Decrease  in  biomass  of  vascu- 
lar plants  is  evident  in  the  second  year  after  introduction  (Figure  1+2), 
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Figure  40.  Simulated  changes  in  total  number  of 
white  amur  for  10  years  after  their  introduction 

T BIOMASS  OF  WHITE  AMUR 


Figure  41.  Simulated  changes  in  total  biomass  of 
white  amur  for  10  years  after  their  introduction 
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Figure  42.  Simulated  changes  in  biomass  of  macro- 
phytes and  epiphytic  algae  for  10  years  after  the 
introduction  of  white  amur 


and  only  10  percent  of  the  initial  amount  of  biomass  remains  by  the 
ninth  year.  This  decrease  in  biomass  slows  the  pumping  of  nutrients 
from  the  sediment,  causing  a decrease  in  orthophosphate  levels  in  the 
epilimnion  (Figure  43),  which  in  turn  leads  to  a decrease  in  phytoplank- 
ton levels  (Figure  44).  Increased  light  availability  permits  the 
epipelic  algae  to  increase  in  biomass  (Figure  45 ) , although  phosphorus 
levels  in  the  hypolimnion  (Figure  46)  are  somewhat  decreased  by  the 
lower  input  from  the  epilimnion  during  turnover. 

91.  Zooplankton  are  much  reduced  by  the  fifth  year  (Figure  47), 
presumably  as  a result  of  the  decrease  in  their  primary  food  source, 
phytoplankton.  Benthic  invertebrates  increase  (Figure  48),  however,  in 
response  to  the  increased  biomass  of  epipelic  algae.  Fish  populations 
as  a whole  increase  toward  the  end  of  the  10-year  period  (Figure  49). 
Herbivorous  fish  do  not  show  a dramatic  change  in  biomass  (Figure  50 ). 
Primary  predator  fish  decrease  somewhat  (Figure  51).  Secondary  predator 
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Figure  45.  Simulated  changes  in  the  biomass  of 
epipelic  algae  for  10  years  after  the  introduc- 
tion of  white  amur 
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Figure  46.  Simulated  changes  in  the  level  of 
orthophosphate  in  the  hypolimnion  for  10  years 
after  the  introduction  of  white  amur 
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Figure  51 • Simulated  changes  in  biomass  of  the 
primary  predator  fish  population  for  10  years 
after  the  introduction  of  white  amur 

fish,  which  feed  heavily  on  benthic  invertebrates  and  other  fish,  in- 
crease significantly  (Figure  52)  as  the  cover  is  removed.  After  an 
initial  rise,  the  amount  of  detritus  decreases  (Figure  53). 
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PART  V : CONCLUSIONS 


92.  The  Lake  Conway  ecosystem  model  as  presently  formulated  con- 
tains those  state  variables  which  are  believed  to  be  essential  for  de- 
scribing seasonal  and  long-term  changes  in  a southeastern  lake  ecosys- 
tem. In  the  final  model,  the  magnitude  of  change  of  the  state  variables 
will  be  determined  by  conditions  in  Lake  Conway,  but  the  nature  and 
rates  of  interactions  among  these  state  variables  should  be  common  to 
lake  ecosystems  throughout  the  region.  Predictions  could  therefore  be 
made  of  the  effect  of  the  white  amur  on  conditions  in  another  lake, 
given  certain  morphological,  climatological,  and  biological  data  about 
the  lake. 

93.  The  current  inadequacy  of  nutrient  relationships  within  the 
model  prevents  the  model  as  it  is  now  formulated  from  being  used  for 
sound  prediction.  Many  of  the  discrepancies  between  data  and  simula- 
tion can  be  attributed  to  treatment  of  relationships  between  nutrient 
availability  and  productivity.  Other  discrepancies,  such  as  the  failure 
of  the  model  to  produce  the  standing  crop  of  macrophytic  plants  that  is 
known  to  exist,  can  be  investigated  experimentally  by  measuring  differ- 
ent growth  parameters.  Correction  of  the  difficulty  with  the  zooplank- 
ton component  can  be  made  by  testing  different  assumptions  in  the  model. 
Other  discrepancies  may  appear  or  be  resolved  when  all  the  data  from  the 
other  research  teams  have  been  provided. 

9b.  It  is  interesting  to  note  the  effects  of  the  white  amur  on 

the  current  ecosystem  model  and  to  speculate  on  the  significance  of  the 

simulated  outcome.  The  increase  in  epipelic  algae  may  be  a phenomenon 

that  is  specific  to  stratified  lakes.  In  a shallower  lake  or  pond, 

epipelic  algae  might  compete  directly  with  phytoplankton  for  nutrients, 

although  they  would  still  be  shaded  by  both  phytoplankton  and  macro- 

88 

phytic  plants.  Singh  et  al.  mention  the  disappearance  of  Spirogyra 
sp.  from  a O.OU-ha  nursery  pond  after  introduction  of  the  white  amur. 

It  is  doubtful  that  the  algae  were  consumed.  Stevenson  mentions  that 
white  amur  did  not  eat  the  Spirogyra  sp.  that  was  fed  to  the  fish. 
Persistence  of  the  algae  in  this  case  was  presumably  due  to  the  continued 
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In  such  situa- 


12 

feeding  of  fish  meal,  which  would  have  sustained  it. 
tions,  epipelic  algae  and  phytoplankton  may  well  he  using  the  same 
resources,  and  neither  will  survive  when  a major  nutrient  source  such  as 
macrophytes  is  removed. 

95.  The  decreases  in  phytoplankton  and  zooplankton  predicted  by 

90 

the  model  are  corroborated  by  Terrell,  who  reported  a decrease  in 
phytoplankton  and  zooplankton  in  acidic , soft-water  ponds  to  which  white 
amur  had  been  added  for  weed  control  tests.  He  attributed  this  to  pre- 
cipitation of  orthophosphate,  iron,  and  magnesium  by  organic  acids. 

96.  Stott  et  al.^1  report  that  bream  in  ponds  stocked  with  white 

amur  showed  faster  growth  rates  than  in  control  ponds.  Production  of 

bass,  crappie,  and  bluegill  in  an  Arkansas  lake  seemed  to  increase  after 

92 

introduction  of  white  amur;  results  for  gizzard  shad  were  ambiguous. 

In  the  same  study,  it  was  noted  that  phytoplankton  and  zooplankton 

blooms  frequently  result  in  eutrophic  situations  when  a high  stocking 

rate  of  white  amur  causes  weeds  to  be  removed  rapidly.  In  cases  such  as 

this,  input  of  nutrients  from  runoff  or  groundwater  seepage  may  be  large 

enough  to  dominate  the  lake  after  the  macrophytes  are  removed  and  no 

longer  exert  shading  or  allelopathic  effects  on  the  phytoplankton. 

Stott  et  al.^1  notes  that  such  a bloom  may  result  "in  extremely  rich 

lakes"  even  when  low  stocking  rates  are  used.  The  level  of  stocking 

itself  probably  has  little  effect. 

93 

97.  Sutton  suggests  that  white  amur  release  the  nutrients  tied 
up  in  aquatic  macrophytes,  making  them  available  for  other  organisms. 
However,  it  appears  more  likely  that  the  fish  remove  these  nutrients 
almost  permanently  from  use.  If  sediment  respiration  exceeds  production, 
long-term  use  of  the  white  amur  may  reduce  the  quantity  of  sediment  on 
which  a crop  of  submersed  macrophytes  depends,  reversing  the  succes- 
sional  trends  that  allowed  the  weeds  to  proliferate  in  the  first  place. 

If  the  assumptions  used  in  the  model  for  nutrient  cycling  are  valid — 
and  many  of  them  are  supported  in  the  literature — the  model  predicts  the 
disappearance  of  sediment,  a parameter  that  is  difficult  to  measure  in 
situ  and  would  require  years  to  observe.  This  will  be  an  important  com- 
ponent to  observe  in  the  Lake  Conway  model  after  it  has  been  improved. 
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Table 


Table  2 


Relative  Productivity  and  Respiration  Rates  for  Submersed 


Macrophytes  at  Three  Different  Temperatures 


* Net  photosynthesis,  based  on  data  for  Hydrilla,  Potamogeton.  and 

Vallisneria.^-9 

**  Respiration  function:  0.63  * e0,0^  * T 
t Gross  primary  productivity  (NP  + R);  values  have  been  converted  to 
a relative  scale. 


Table  h 

Food  Preferences  of  Adult  Primary  Predator 
Fish  Found  in  Lake  Conway 


,1 

$ 


! 


Fish 

Fundulus  seminolis 
( seminole  killifish) 

Luc an i a goodei 

(bluefin  killifish) 

Labidesthes  sicculus 
(brook  silverside) 


Enneacanthus  gloriosus 
(blue-spotted  sunfish) 

Lepomis  gulosus 

(warmouth  sunfish) 

L.  macrochirus 
(bluegill) 


L.  microlophus 
(redear  sunfish) 

L.  marginatus 

(dollar  sunfish) 

Etheostoma  fusi forme 
(swamp  darter) 

Noturus  gyrinus 
(tadpole  madtom) 

Heterandria  formosa 
(least  killifish) 


Fundulus  chrysotus 
(golden  topminnow) 

Gambusia  affinis 
(mosquitofish) 

L.  punctatus 

(spotted  sunfish) 

Ictalurus  natalis 
(yellow  bullhead) 

Aphredoderus  sayannus 
(pirate  perch) 


Major  Food  Items References 

1*0 

Benthic  invertebrates,  seeds  McLane 
of  aquatic  plants 

1*9 

Benthic  invertebrates  McLane 


Insects  and  crustaceans 


Insect  larvae,  small 
crustaceans 

Benthic  invertebrates,  small 
fish 

Zooplankton,  benthic  inver- 
tebrates, crustaceans, 
aquatic  plants 


Benthic  invertebrates 


Ewers  and 
Boesel,^ 
McLane^-9 


Chable,53 

McLane 


Hubbell-^ 


Chable,53 
Emig,55 
Flemer  and 

Woolcott ,56 

Huish5T 

Emig,58  Huish5' 


Zooplankton,  benthic 
invertebrates 

Benthic  i nvert ebr at e s , 
zooplankton 

Zooplankton,  benthic 
invertebrates 

Benthic  invertebrates , zoo- 
plankton , filamentous 
algae 

Benthic  invertebrates, 
zooplankton 

Insects,  zooplankton,  fila- 
mentous algae 

Insects , crustaceans 


Chable,53 

Me Lane ^9 

Me Lane 

Carlander , 
McLane^9 

McLane^*9 


McLane^ 


McLane**9 

53 

Chable^ 


Benthic  invertebrates,  zoo- 
plankton, aquatic  plants 

Benthic  invert ebrat s , zoo- 
plankton, aquatic  plants 


1*9 

McLane , 
Miller59 

1*9 

McLane 


Table  5 

Food  Preferences  of  Secondary  Predator 
Fish  Found  in  Lake  Conway 


Fish 

Major  Food  Items 

References 

Lepisosteus  platyrhincus 
(Florida  gar) 

Fish,  crustaceans,  insects 

UU 

Carlander , 

Crumpton^O 

Esox  niger 

(chain  pickerel) 

Fish 

„ . 6l 

Buntz , 

Carlander^1* 

Micropterus  salmoides 
(largemouth  bass) 

Fish,  crustaceans,  benthic 
invertebrates 

62  T h9 
Emig , McLane 

Pomoxis  nigromaculatus 
(black  crappie) 

Fish,  benthic  invertebrates 

6 3 

Goodson, 

Huish^T 

Lepisosteus  osseus 
(longnose  gar) 

Fish,  benthic  invertebrates 

_ ^ 6o 

Crumpton , 

Lagler  and 
Hubbs*^ 

Ami a calva 
(bowfin) 

Fish,  benthic  invertebrates, 
zooplankton 

Lagler  and 
Hubbs6{* 

McLane^9 

Esox  americanus 
(redfin  pickerel) 

Benthic  invertebrates,  fish 

Lq 

McLane 

Table  6 

Phosphorus  Budget  for  Lake  Conwa; 


h 

P 


i 


Month 

Rainfall* 

cm 

Total  Phosphorus 

2 

g/m  , from 

Input** 

Total  Monthly 
Input  from  All 
Sources 

2 

g TP/m  per  mo 

Rainfall 

Citrust 

Runoff 

Urban ft 
Runoff 

January 

12 

0.0053 

0.0023 

0.0297 

0.0373 

February 

7 

0.0031 

o.ooii* 

0.0175 

0.0220 

March 

11 

0.001*8 

0.0021 

0.0273 

0.0342 

April 

7 

0.0031 

0.0011* 

0.0175 

0.0220 

May 

12 

0.0053 

0.0023 

0.0297 

0.0373 

June 

17 

0.0075 

0.0033 

0.01+23 

0.0531 

July 

16 

0.0070 

0.0031 

0.0395 

0.0496 

August 

18 

0.0079 

0.0035 

0.01*1*7 

0.0561 

September 

29 

0.0128 

0.0056 

0.0720 

0.0904 

October 

3 

0.0013 

0.0006 

0.0073 

0.0092 

November 

2 

0.0009 

0.0004 

0.0049 

0.0062 

December 

7 

0.0031 

0.0011* 

0. 0175 

0.0220 

TOTAL 

l4l 

0.0621 

0.0271+ 

0.3499 

0.4394 

* Precipitation  data  from  Reference  73. 

**  Shannon  and  Brezonik  report  an  average  concentration  of  0.044  g 
total  phosphorus/m3  in  rain."^  For  runoff  from  citrus  and  urban 
lands,  they  report  average  total  phosphorus  loadings  of  0.018  and 
0.11  g/m2  per  year,  respectively.  Monthly  loading  rates  from  land 
sources  were  determined  by  using  the  following  formula: 

Monthly  increment  = (area  of  land  type  in  basin) 

x (yearly  loading  rate)  x (percent  of  total  pre- 
cipitation which  occurred  during  that  month) 

This  equation  assumes  that  runoff  from  the  land  is  immediate  and 
directly  proportional  to  the  rain  which  falls  upon  it. 
t The  area  of  citrus  lands  in  the  drainage  basin  is  11. 254  x 10°  m2 
(E.  Blancher,  personal  communication), 
ft  The  area  of  urban  lands  in  the  drainage  basin  is  23.47  x 10°  m2 
(E.  Blancher,  personnal  communication). 


Table  7 

Growth  Efficiencies  of  White  Amur  of  Different  Ages  and  Weights 


Growth 

Food  Intake* 
percent 


80-90 

100-199 

200-299 

b!9 

680 

700-799 

955 


3.h-k. 3++ 
1.2-2. 8 
2. 6-3. 8 
0.7 
0.1 

1. 0-2.1 
1.0 


* Percentages  are  calculated  from  wet  weights  of  both  food  and  fish. 
**  Data  summarized  in  Fischer  and  Lyukhnovich. 9 
+ Data  from  Singh  et  al. ° 

ft  Ages  estimated  from  data  summarized  by  Fischer  and  Lyakhnovich. 


Table  8 

Daily  Weight-Specific  Digestion  Rates  of  White  Amur 


Initial  Wet 
Weight,  g 


Daily  Food 
Intake , g 

115-1350 

500-2300 


Intake 

Biomass,  percent 


12-135* 

25-115* 


30-130** 


* From  data  summarized  by  Fischer  and  Lyakhnovich. ^ 
**  Stott. 89 
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APPENDIX  A:  PROGRAM  LISTING  FOR  THE  LAKE  CONWAY  MODEL 


INITIAL 

♦STATE  VARIABLES 
INCON  ICFN=9.2lE-4 
INCON  1 CEB—  *4182 
INCON  I CB 1=1*73 
INCGN  ICDT=150. 

INCON  ICE  P= I • 

INCON  ICHF=.07 
INCON  ICPEPI=.006 
INCON  ICPF1=.74 
INCON  ICPF2=.36 
INCON  I CPHYP= «0  05 
INCON  ICPP=.84 
INCON  ICTJ=30. 

INCON  ICVP=45. 

INCON  ICZP=.21 

♦BENTHIC  INVERTE8RATE  CONSTANTS 
CONST  KBI 1 =1 • 

CONST  KBI2=8.28E-3 
CONST  KB 13= 1 • 26 
CONST  KBI4=.4 
CONST  KBI 5= • 1 
PARAM  bIASSM=.689 

♦UETRITUS  CONSTANTS 
CONST  K J=2 1 • 75 
CONST  KK=10.87 

♦EPIPELIC  ALGAE  CONSTANTS 

CONST  K IE P=  31  3*  48 

CONST  K2EP=3.2 

CONST  K5EP=1.34 

PARAM  EPMEFF=. I 0 

PARAM  EP50=350.0 

♦ HE  RBI VOROOS  FISH  CONSTANTS 
CONST  KHF1=.462 

CONST  KHF2=.183 
CONST  KHF3=3.41E-J 
CONST  KHF4=.153 
CONST  KHF5=.121 
CONST  KHF6=1. 

CONST  KHF  7= 1 .5 
CONST  KHF8=.l 
CONST  KHF9=.I 
PARAM  HFASSM=.69 

KPHHF 1=2 .44E-4 

♦HYPOLIMNET IC  PHCSPHORJS  CONSTANTS 
CONST  KPHYP1=.06 
CONST  KPHYP2=100* 

PARAM  PEPX=.009 
PARAM  PHYPX=.005 
CONST  KC2P=.01 
CONST  PCTPIF=.S 

♦KC2PAPCTPIF  IS  THE  AMOUNT  OF  NON  REFRACTORY  P IN  A FIaHJKC2P  CONVERTS 

♦ FISH  CARBON  TO  PHOSPHORUS;  50  PERCENT  OF  THIS  IS  AVAILABLE  FOR 

♦ REMINERALIZATION 

♦LIGHT  FUNCTION  CONSTANTS 

CONST  KC=.949 

CUNST  KD=-.456 

CONST  KI=-.28 

CUNST  KM=-.43 

PARAM  REFLCT=.05 

CONST  Z=2.0 

♦PHYTOPLANKTON  CONSTANTS 
CONST  K1PP=372. 

CONST  K2PP=.76 
CONST  K4PP=36. 


I 


M 

I 

CONST  K6PP=10. 

CONST  K7PP=.  1 95 
CONST  KBPP=.a05 
CONST  K9PP=C. 

PARAM  PPMEFF  = . 0098 
P ARAM  PP50=14597» 

* THE  AMOUNT  OF  PHOSPHORUS  RELEASED  FROM  DE  AO  PHY  TOPL  ANnT  ON  IS  EOUAL  TO 

* 75  PERCENT  OF  THAT  AVAILABLE.  THE  OTHER  25  PERCENT  IS  ACCOUNTED  r OR 

* IN  ZOOPLANKTON  EXCRETION 

♦PRIMARY  PREDATOR  FISH  CONSTANTS 
CONST  KPF 1 1 = • 85 
CONST  KPF  12=7. 07E-3 
CONST  KPF 1 3=4. 

CONST  KPF 1 4= .318 
CONST  KPFI5=1.7E-2 
! CONST  KPF16=3.8 

CONST  KPF I 7=  » 1 
PARAM  P 1 A SSM- • 8 
CONST  KPHP 1 1=2. 44E- A 

♦SECONDARY  PREDATOR  FISH  CONSTANTS 

CONST  KPF21 = .298 

CONST  KPF22= • 166 

CONST  KPF23=2.98E-2 

CONST  KPF24=.l 

CONST  KPF25= .6 

CONST  KPF  26=  5. 95E— 2 

PARAM  PF2ASM= 1 . 

♦TUBER  CONSTANTS 
CONST  KTUl=l .5 
CONST  KTU2  = 1 6 

♦VASCULAR  PLANT  CONSTANTS 
CONST  KlVP=2.2E-2 
CONST  K3VP=. 58 

| CONST  K6VP=.I 

CONST  K7VP=.05 
CONST  K8VP=. 65 
CONST  K9VP=.7 
CONST  KCAL2C=. 1 1 
PARAM  VPMEFF=. 1426 
PARAM  VP50=2726 .48 

♦♦HITE  AMUR  CONSTANTS 
CONST  K1 FN= .0833 
CONST  K 1FB=I . 94 
CONST  K2FB=1 .0 
CONST  K3FB= .6065 
CONST  KW2C=. 036 
CONST  KPHFB 1=2 .44E-4 

♦YOUNG  PRIMARY  PREOATOR  FISH  CONSTANTS 
CONST  KYNG 1=3. 53E-2 
CONST  KYNG2=1.4E-2 
CONST  KYNG 3=2.62E  — 4 
CONST  KYNG4=3.37E-2 
CONST  KYNG 5 =8 . 96E— 3 
CONST  K Y NG  6 = • 1 06 
CONST  KYNG7=8. 

CONST  KYNG8=3. 76E-2 
CONST  K YNG9=.  1 
PARAM  YNGASM= .69 

KPHY  G I =2 . 44 E — 4 
♦ZOOPLANKTON  CONSTANTS 
CONST  K ZP l = 4 • 

CONST  KZP2=9.57 
CONST  KZP4=. 142 
CONST  K ZP5=6 • 5 
PARAM  ZPASSM=.47 


♦OTHER  CONSTANTS  AND  PARAMETERS 
PARAM  KO XI C =6 . 4 
PARAM  KOVR 1=3.01 
PARAM  KOVR=l .37 
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.1. 

.4.2 

• 2*  »4  *3i 3* 

,4.4 

.4 

4.5.5  .-#4 . 6 

.6. 

.4  .7 

• 7 • • 4 «8  *8  • 

.4.9 

.9 

5.0.0..S. 1 

. 1. 

.5.2 

•2* • 5 • 3 • 3 • 

.5.4 

.4 

5.5. 5 • . 5 . 6 

• 6. 

,5.7 

» 7*»5«8»8# 

. 5.9 

. 9 

6.0.0. .6 . 1 

.1. 

.6  .2 

• 2*  »6 *3*3* 

.6.4 

. 4 

6. 5.5. .6. 6 

.6. 

.6.7 

• 7 • • o • 0 t 0 • 

.6.  9 

.9 

7. 0.0. -.7.1 

.1. 

.7.2 

• 2 • • 7 • 3 • 3 • 

.7.4 

. 4 

7.5.5. *7.0 

• 6. 

.7.7 

•7* *7*0*0# 

.7.9 

.9 

0. 0.0. .3.1 

. 1. 

.8.2 

•2**8*3*3* 

.8.4 

. 4 

8. 5. 5. .0.0 

, 6 . 

.8  .7 

•7**8*8*0* 

.8.9 

. 9 

9.0.0.. 9. 1 

. 1. 

.9.2 

• 2*  *9*3 *3* 

.9.4 

.4 

9.5. 5.. 9.0 

.6. 

.9.7 

•7**9*8*b* 

.9.9 

.9 

1 0.0.0-.1 0 

.1  . 

1 . 

FUNCT  SN  = 0 . .61 . 4E3.  .85.  71.40E3.1.64.134.8E3 


4.16.1 13.CE3 .3.0. 150.0E3. 5.34. 143.0E3.6.69. 113. 9E3. 7.51 
108. 0E  3.8. 33. 71. 4E 3. 9.  15.93.GE3.10..61.4E3 


FUNCT  DG 1=0. ..58. .53. .55. 1.38.  .85. 2. 3, .84, 3 

5. 74. .  80  0. b. 39. .950. 7. 38. . 97. 6. 15 

FUNC  T DG2=0. . .20.. S3 ..2  0 . 1.38. .25. 2. 3. .2 5. 3 

5.74. . 670 .6.3  9. .5  70.7.  38..  50 . 8.  15 
FUNCT  DG3=0. . .92. .53. .90. 1 . 38 . 1 . 0 . 2. 3. . 99 . 3 

5. 74. . 92 0.6. 39. .950. 7. 38. .96. 8. 15 

FUNCT  DG4=0. • .64. .53. .52. 1.38. .76. 2. 3. .78. 3 

5. 74. 1.06. 6. 39. 1.03.7.33..  98. 8. 15 

FUNCT  DG6=0.  ..36.. 5 3. .34. 1.38. .52. 2. 3. .55. 3 

5.74.1.24.6.39. 1.03.7 .38. . 97. 6. 15 

FUNCT  ORTHO—  0 .. .0145.. 85. .00 86.  1.64,  .0133.2 
.0207,5.0.. 0193.5.84 ,.0219 .b.69 
. 00  24 . 9. 15.  .0  086 .10.0. .0145 


2. 5. 165.  05 J.3.32.84.00E3,... 


04,1  .0.4.07,  .99,  4. 95. .90,... 
.88, 9.93, .65, 10.  ,.56 

04..  47. 4 . 07.  .49, 4 . 95, .60  .... 
.32,9.93,. 21  .10.  , . 20 

04. . 98.4. 07.. 97. 4. 95. .95.... 
.99, 9. 93. .94. 10. ..92 

04. . 91. 4. 07. .94.4.95.1.0.... 
•80, 9. 93,. 65, 10. « .64 

04. . 82.4. 07. .87.4. 95. 1 .1  .... 
.60. 9.93. .38 , 10 . , .36 

5,  . 0036.  3.32.. 0145.4. 16... . 
.0353.7.51. .0036. 8.33.. • • 


FUNC  T PE  RCNTs  0. 0 ». 85 ,1000. ,.888, 2000. ,.915, 3000. ,.94.4300. ,.96,5000.,... 
.97. 6000. ..98. 7000. ..99. 80 00... 995,9000. ,.99 8, 1.54,1 . ,1 .E3.1. 

F UN  CT  2P  T — 1 0 • ,.651 ,1  .53 ..651 . 11. 33, .771 , 12. 3, .76  1.  13.04.  .891  , .. . 

14. 07. *911. 14. 95, .971, 15. 74, 1.1 ,16.39, .941 ,17. 38. .921.... 

18. 15.  *791. 19. 93. .67 1.110. ,.671 

FUNCT  BIT— 10...451.1 .53 ,.451.(1.38, .581, 12.3. .71, (3.04, .851 •••• 

14.  07.  *891, 14. 95. .96 1,15. 74. 1.1.  16. 39.  .98). 17. 38, .96}.... 

18. 15.  *771. 19. 93. .51 1.110.  ,.51  1 

FUNCT  FSHT=10. ..35}.t.53..351.tl.38..67}.t2.3..7}.l3.0«,.91}.... 

14. 0 7..931. 14. 95.. 98 1.15. 74. 1.1. 16. 39, .971, 17. 38,  .96  1,... 

18. 15,. 75 1.19. 93.  .411, 110.  ..41) 

FUNCT  HEIGHT  = 1 0. 0.2. 01. 15. 18. 1.91, 19. 69.1. 81. 112. 7 8, 1.71. 115. 3.1. bl.... 

117. 6 7. 1. 51. 1 19. 8 2. 1. 41. 121. 84. 1.31. 123. 74. 1.21.125. 31. 1.11.. .. 

1 26.52.1 . 1. 127 .73. .  9 1 . t 28 . 52 ,. 8 1 , 1 28 .97,  . 7 1 . 1 29 . 32 .. 6 1 ... . 

129. 51. . 51. 129. 58. .41. 129. 6. .31. 12 00. ..2) 

FUNCT  SHAOE=(0..1 .1.1150., .01 1 
FUNCT  BLOCKS  10 . . 1 . 1 . I 150 . , .51 

FUNCTION  PLANTS=10.  .2.  1 . 1 45 . . 1 . 1 . 1 1 3 0 . . 1 . 1 . 1 20 0 . . .5 1 

FUNCT  ORGCC1  = 10 . ..0 1 1.115..  .141.  130. ,.181. 145., .141,1  60. ..Oil. 1100. ..Oil 
FUNCT  ORGLC2sl 0.,.011.t5.,.14!,tlO.,.18!.tl6...14|,(20...Ql).tlGG.,.Ull 


DYNAMIC 


T I MEX  = AFGEN1 DC V. T I ME  1 

SUN=AFGEN1SN.T 1MEX1 

OPI  NPTsAFGENl  CRT-tO.T  I MEXl 


A3 


T1PP=AFGEMDG1  .TIME  X) 
T2PP=AFGEN(uo2. T J MEX) 
TIEP=AFGEN(0G1 .T IMEX) 
T 2EP=AF GE  N<OG 1 . TI  HEX) 
TIVP=AFGEN(DG3.T  IMEX) 
T 2 VP— AF  GE  NT  0G4, TIMEX) 
T2NC=AFGEN(DG0. TI  ME  X) 


ZPRI ME=AFGEN( HE IGHT  , VP ) 

RATI 0=AFGEN( PERCY  T , I ) 

VPSH01=AFGEN( SHADE. VP ) 

VPSHAD=AF GEN (BLOCK. VP) 

2PT  EMP= AFGEN ( ZPT » T I M EX ) 

FSHTMP=AFGEN(FSHT, TIMEX) 

BIT  EMP=  AFGENI B 1 T»  T IMEX) 

CPPRLS=AFGEN( QRGLC1 .GPPPP) 

CEPRLS=AF GE  N ( 0RCLC2 , GPPEP ) 

* THE  PLANT  EQUATIONS  INCLUDE  THE  FOLLOWING  VARIABLES 

* GPPXX  IS  GROSS  PRIMARY  PRODUCTIVITY 

* XXRESP  IS  RESPIRATION 

* XXHERB  IS  LOSS  OF  BIOMASS  DUE  TO  HERB  I V OR  Y 

* XXLEAC  IS  LEACHING  OF  DISSOLVED  ORGANIC  CARBON 

* CXXRLS  CONTROLS  ACTIVE  SECRETION  OF  ORGANIC  COMPOUNDS.  IT 

* IS  GREATEST  AT  INTERMEDIATE  LEVELS  OF  PRODUCTIVUY. 

♦SUBMERSED  VASCULAR  PLANT-EPIPHYTIC  ALGAE  COMPLEX  (VP) 

VP=INTGRLt ICVP.GPPVPFSPROUT -V PRESP-SLOUGH-PRFCRV -VPHER d ) 

GPPVP  = K 1 VP*VPLI T£  *VP*(  VPMEFF/I  1.FVPLITE/VP50))*T1VP*KCAL2C 
VPL I TE— SUN* ( 1 . — REFLCT ) * < XC« EXP < KD» ZPR IME FK  I*ZPP  IME  *PP/2  ) ) 

*VPL ITE  IS  THE  LIGHT  AVAILABLE  FOR  VP  PHOTOSYNTHESIS  AnD  IS  A 

* FUNCTION  OF  THE  HEIGHT  OF  THE  VP  AND  THE  OPTICAL  PROPERTIES  OF  N20 
•VPMAX  IS  THE  MAXIMUM  UBSERVEO  UIOMASS  OF  vP 

♦ VPMEFF  IS  THE  MAXIMUM  PHOTO SYN THET I C EFFICIENCY  OF  VP 

* VP50  IS  THE  LIGHT  INTENSITY  AT  WHICH  SO  PERCENT  OF  VPMEFF  OCCURS 
•REFLECT  IS  THE  DECIMAL  PERCENT  OF  SOLAR  ENERGY  WHICH  IS  REFLECTED 

* AND  CHANGES  SEASONALLY 

* ZPR I ME  CORRELATES  CUMULATIVE  VP  BIOMASS  WITH  THE  DISTANCE  FROM  ThE 

* TOP  OF  THE  CANOPY  TO  THE  WATER  SURFACE)  THE  SECOND  TERM  IN  THE 
» EXPONENT  CALCULATES  SHADING  OF  VP  BY  PHYTOPLANKTON 

VPRE  SP=K3 VP*  VP*T2 VP 

VPHERB=BIXVPFBFXVP+YNGXVPfCPXVPFKw2C*FCNSMP 

* SLOUGH= SLOUGH I NG  OF  LEAVES  AND  OTHER  PLANT  PARTS 
SLOUGH=  K6  VP*VPFK9 VP*VP*OIE 
FALC=IMPULS (.75.1.) 

D1E= PULSE (.25 .FALL) 

•P04-P  LEACHING  BY  VP  (VPPLCH) 

VPPLCH=KC2P*VPRESP* .5 

•SLOUGHED  MATERIAL  WHICH  REACHES  THE  SEDIMENTS  (VP2SE3) 

VP2SED=Ke VP* SLOUGH 

•PU4-P  RELEASED  UPON  DECAY  OF  SLOUGHED  MATERIAL  IN  H20  COLUMN  ( PSLOF ) 
PSLUF=KC2P*  K7 VP *S LOUGH 

* FORMAT  TON  AND  GERM  INAT ION  OF  PROPAGOLES  AFFECT 

* VASCULAR  PLANT  BI GMASS  IN  THE  SPRING  ANO  FALL 

•TUBERS  (TUI 

TU= I NT GRL( TCTD.PRFORM-SPROUT ) 

PRFORM=DIE*KTUl»VP 
SPROUT=BIRTH#KTU2  *TU 

•PHYTOPLANKTON 

PP=I NTGRL (ICPP .GPPPP-P  PRtSP-PPLt AC-PPHERB- SINKNG) 

GPPPP=K1 PPWPPLI TE* (PPMEFF/ U .fPPL 1TE/PP5  0 ) ) *Pt  P I » T 1 PP  *K  CAL  2C*PP 

•PPLITE  IS  THE  LIGHT  AVAILABLE  FOR  PP  PHOTOSYN THESI S ANO  IS  A FUNCTION 

* OF  THE  OPTICAL  PROPERTIES  OF  H20,  SELFSHADING.  AND  iHADlNG  BY  VP 
•PPMEFF  IS  THE  MAXIMUM  PHOTOSYNTHETIC  EFFICIENCY  OF  PHYTOPLANKTON 
•PP50  IS  THE  LIGHT  INTENSITY  AT  WHICH  SO  PERCENT  OF  PPMEFF  OCCURS 
•THE  EXPOFfcNT  TERM  REPRESENTS  SELF  SHADING! VPSHAO-SH ADI NG  BY  VP 

PPL  I TE=SUN*( 1 .-REFLCT ) »K2PP*txP( KM *P P) * VP SHAD 


Ah 


I 

\ 

I 

PPRE  SP=K4PP*T2PP*PP 
PPLE AC=CPPW_b*&P=PP 
PPHERb=ZPXPP*HFXPP+YNGXPP 

•SI NKNG=SETTL ING  OUT  Or  PHYTOPLANKTON 
SI NKNG=K6PP*PP 

*PP2SED=THE  PERCENTAGE  OF  SINKING  PH YT UPLA NK TON  WHICH  REACHES 
•THE  SEDIMENT 

PP2SED  = K6PP*  SI NKNG 

•P04-P  RELEASED  UPON  DECAY  OF  SINKING  Ph YTOPL A NK TON= PS 1 NK 
PSI NK=KC2P* (SI NKNG-PP2SLD I *.5 
PPPART=K7PR*Kt>PP*PP 


♦EPIPELIC  ALGAE  (EP) 

EP= I NT  GRL < ICEP.&PPEP-t PRE SP-EPLE AC-EPHERB-EPMORT ) 

GPPEP=KlEP»EPLITE»l EPWEFF/ ( 1 . ♦ EPL I TE/E PS 0 > > »PH  YP^*  T 1 F T> » KC AL2 C* LP 

•EPMEFF  IS  THE  MAXIMUM  PHOT OSY NT HET I C EFFICIENCY  OF  EP 

•EP50  IS  THE  LIGHT  INTENSITY  AT  WHICH  SO  PERCENT  OF  PPWEFF  OCCURS 


• EPLITE  IS  THE  LIGHT  AVAILABLE  FOR  EP  PHOTOSYNTHESIS  A V D IS  A 

• FUNCTION  OF  THE  OPTICAL  PROPERTIES  OF  H20.  AND  SHADING  BY  VP  ANd  PP 

EPLI TE=SUN* ( 1 • — REFLCT ) *EXP  <KD* ZFKM  WPP ) *VPS HD  I 

EPRESP=K2EP*T2EP»EP 
EPLEAC=CEPRLS*GPPEP 
EPHERB=ZPXEPfHFXEP+PF1  XEPf  yngxep 

• EPMORT=MORT AL I TY  OF  EPIPELIC  ALGAE 

EPMORT=K5EP*EP 

•THE  ANIMAL  EQUATIONS  INCLUDE  THE  FOLLOWING  VARIABLES 

• XXA  SSM*  XXF  OOO  (OR  XXX ASM*  XXXF D J IS  THE  AMOUNT  OF  FOOD  ASSIMILATED 

• XXRESP  (OR  XXXRSP)  IS  THE  METABOLIC  RATE 

• XXPRED  (OR  XXXPRD)  IS  THE  LOSS  DUE  TO  PREDATION 


•THE  ABILITY  OF  A PREDATOR  TO  CAPTURE  ITS  PREY  IS  IN  PART  DEPENDENT  ON 
• THEIR  RELATIVE  ABILITIES  TO  FUNCTION  AT  THE  AMBIENT  TEMPERATUkE 
BI ZP— BI  TE  MP/ZP TEMP 
FSHB I=FSHTMP/BI TEMP 
FSHZP=FSHTMP/2PTEMP 


•BENTHIC  INVERTEBRATES  (BI) 

B I = I NTGRL ( I CBI  .BI A SSM* d IF 000 -fa  I RESP-B I PR EO ) 

B IFQQD— B IXZPFBIXVPFBIXOT 
B I E A TN=8 1 ASSMWBI FOOD 
BIXZP=KBIl*ai*ZP*6IZP*BITEMP 
BI XVP=KBI2*BI * VP* 6 I TEMP 
BIXDT=K8I3*BI *B I TE  MP 

BIRESP=KBI4*BI*BITEMP 
BI  PREO=  HF  XBI  + PF  1XB  I *Y  NGXB I -FPF2  XBI 
BIEGEX=(1 .-BIASSM )*BIFOOO 
PHBIEE=KC2P*BIEGEX 

•HERBIVOROUS  FISH  IHF) 

HF=INTGRLI ICHF .HF ASSM*HFFUUO-HFRES P-HFPREO ) 

hffood=hfxepfhfxpp*-hfxvpfhfxdt  fhfxb  i fhf  xzp 

HFEATN=HFASSM*HFFOUO 
HFXEP=KHF 1*EP*HF*FSHTMP 

HFXPP=KHF2*PP*HF*FSHTMP 

HFXVP=KHF3#VP*HF*FSHTMP 

HF  XO  T = KHF  4WHF  *F  SH  TMP 

HFXB  I = KHF5*HF  *BI  *FSHb  I *FSHTMP 

HFXZP=KHF6*HF  *ZP*FSHZP*FSHTMP 

HFRESP  = KHF7*HF  *FSHT  MP 
HFPREO=KPF22*HF*PF2*FSHTMP*COVER 
HFEGEX  = ( 1 . — HF  A SSM ) * HF  ODD 
PHHFEE=KPHHF 1 *HFEGEX 


A5 


JJ 
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♦PRIMARY  PREDATOR  FISH  ( PF 1 ) 

PF1-1  NTGRi-I  I CPF  I . PI  ASSMaPF  1 FDF GKO* TH-bREED-PF  1 PR U-PF 1 R SP  > 

PF IFD=PF l XBIfPF 1 X VPFPF 1 XZP *PF 1 XEP 
PF 1 E TN=P1 ASS  M A PF I FO 

PF1 XHI=KPF I 1 Afll A^F 1AFSHBIAF  SHTMP 
PF1XVP=KPF12*VPAPF1 AF  SHTMP 
PF 1 XZP=KPF1 JaZPAPFI AFSHTMPAFSHZP 
PF l XEP=KPF 1 4 AEPaPF 1 AF  SHTMP 

PFIRSP=KPF1 6APF1 AF SHTMP 
PF 1PR0=KPF21 APF2APF 1 AF SHTMP AC OVcR 

PF1EGX=<1 .-PIASSMJAPF1FD 
PHPlEE=KC2RAPFlfcSX 

BREED  AND  GROWTH  REFER  TO  REPRODUCTION  AND  RECRUlTMt_NT  OF  PRIMARY 
PREDATOR  FISH  . 

GROWTH=K YNG7  AY NG*F SHTMP 
BREED=XPF I 5 APF 1 AB I RTH/DELT 
BIRTH=PULSE( . 1 .SPRING! 

SPR1 NG= I M PULS  1.3.1.) 

♦YOUNG  PRIMARY  PREDATOR  FISH 

YNG=INTGRL<  0 . . BREED FYNGASM AYNGF D-GRO A TH- YNGRSP— Y VGPRD ) 

yngfd=yngxepfyngxppfyngxvpfyngxdtfyngxbi  fy  ngxzp 

YNGSTN= YNGASMA YNGFD 
YNGXEP=KYNG1 AEPAYNGAF  SHTMP 
YNGXPP=K YNG2APPAYNGAF  SHTMP 
YNGX VP=KYNGJA V PAY NGAF SHTMP 
YNGXDT  = KYNGAA  YNG»F  SHTMP 
YNGXBI=KYNGSAYNGABIAFSHTMPAFSHB I 
YNGXZP=KYNG6A YNGA ZPAF SHTMP AFSHZP 

YNGRSP=KYNG8AYNGAFSHTMP 

YNGPRD=KPF23aPF2A YNGA F SHTMP *Cu VER 
YNGEGX-=  ( 1 .-YNGASM ) A YNGFD 
PH  YNGE  = KC2P*  YNGEGX 
ALLHF=PF1 FYNG 
ALLF=HFFPFIFPF2FYNG 

♦ SECONDARY  PREDATOR  FISH  < PF  2 ) 

PF2=I NTGRL< I CPF 2, PF2 AS MaPF 2FD-PF 2RSP-PF2DT  H) 

PF2FD=PF 1PRDFHFPREDFYNGPRDFPF2XB1 
PF2XBI=KPF24aBI APF2 AFSHTMPAFSHBI ACOVER 
PF 2ETN=PF  2ASMA PF2FD 
COV  ER=  AFGENI PLANTS. VP ) 

PF2RSP=KPF25APF2AFSHTMP 

♦ PF2DTH  IS  THE  DEATH  RATE  OF  SECONDARY  PREDATOR  FISH 

PF  2DT  H=KPF2to  APF2AFS  HT  MP 

♦ZOOPLANKTON  IZP) 

ZP= I NT GRL IICZP .ZRASSMAZPFOOD-ZPRESP-ZPPRED ) 

ZPFOOD=  ZPXEPFZPXPPFZPX VP 
ZPEATN=ZPASSMAZPF ood 
ZPXEP=KZP1 AEPAZPAZPTEMP 
ZP  XPP= KZP2APPA  ZP ♦ ZP  TEMP 
ZPXVP=KZP4*VPAZPAZPTEMP 

ZPRESP=KZP5AZPAZRTEMP 
ZPPRE  D=  B I XZPfPFIX  ZPfhFXZPf YNGXZP 
ZPEGE  X=  ( l.-ZPASSMIAZPFOUD 
PHZPEE=KC2PAZPEGEXA ,2b 

♦EPILIMNETIC  PHOSPHORUS  EQUATIONS 

PEP I = INT GRL ( I CPEPI • OP  I VPT F V PPlCHFP SLOP FKC2P ♦ I PPL t AC F PPRE SP» . 5-GPPPP ) . . 

fphzpeefphhfeefphpi  EEFPHYNGE FPHF6EE FTRNOV-«a(  PHYPO— PEP  1 ) ) 

♦HYPOLIMNET IC  PHOSPHORUS  EQUATIuNS 

PHYPO=I NTGRL ( I CPHY  P . KC2  P ♦ KK ♦ T 2 NC AANOX I CFKPHY P 1 AK C2P AK J aT 2NC ♦ AEROUC. . . 
FPSLUF  FPSI NK-KC2PAI GPPEP- (EPRESPA .b I-EPLEAC- . . . 

< EPMORTA .3 ) ) FpHblEE  FKC2PAPCTP IF ♦( PF 2DTHFKnF 9AHF RE  SP. . . 
FKPF17APF1  RSPFYNGRSPFFMURT  ) FT  RNOVRA (PEP  I-^HYPU ) ) 


*«HITE  AMUR  EQUATIONS 

*FN=FISH  NUMBER. Fto=TOTAL  FISH  BIOMASS.  l=olOMASS/ INDIVIDUAL, 

*FC NSMP=C QNSUMPTI uN  OF  AQUATIC  PLANTS  BY  THE  AH  I TE  AMJR 
*FASSIM=PERCENT  OF  CONSUMPTION  »HICH  IS  ASSIMILATED 
*FRESP=RESPI RAT  ION  OF  THE  WHITE  AMUR 
*F  MORT  = D I E— OFF  OF  THE  WHITE  AMUR 

FN=INTGRL( ICFN.-K1FNWFN) 

FB=I NTGRL (ICFB.FASS IM-FRESP-FMORT ) 

F CNSMP=K  1FBAA.2FB*  FB*  VP*F  SH  TMP 
PASS IM=K3FB»FCNSMP 
FEGEX=KW2C*(FCNSMP-FASSIM ) 

PHFBEE=KPHFB1 WPEGEX 
FRESP=K5FO#FB*FSHTMP 
FMO«T=Kl F N*FN* I 

•VARIABLE  COEFFICIENT  XSFB  ADJUSTS  RESPIRATLRY  LCSSES  ACCORDING  TO  THE 

* RATIO  OF  ( ( FRESPFFMDRT l/FASS I M ) . THE  RATIO  VARIES  BlT.EEN  0.U5  Fuk 

• A 0-2  YEAR  OLD  FISH  AND  1.0  FOR  A 5-6  YEAR  OLD  FISH 

K5FB=(RAT 10*FASSI M-FMORT 1/ <FB*FSHTmP) 

I =F B/FN 

AN0XIC=ANOlAN,0XICI 
AN=CGMP AR  IT  IMEX. 2 .2  ) 

OXIC=COMPAR(KQXlC.TIMEX) 

AEROBC=NANO< AN, OX IC ) 

TRNOvR=IOR(ISTRMl  .ISTRM2) 

I STR  M 1 = ANU< TKN1 ,0  VR 1 ) 

TRNI=COMPAR( T IMEX , 2.4* ) 

OVRI =COMPAR( KOVR1  .TIMEX) 

I STRM2=EOR<  TRN.OVR ) 

T RN— COMP ART  T IMEX, 7 .5 ) 

O VR— C OMP AR4  KO VR ,T IMEX) 

•DETRITUS  EQUATIONS 

DT  = l NTGRL (I  COT ,. 9*  SL0UGH+PP2SE0+ . 75*ZPEGEX +B IEGEX  + HFEGEXF... 

PF lEGXEYNGEGXFKO I5*BIRE  SPFKHF  8WHFRE  SP  t-KPF  17*P=1SSP*... 
KYNG9WYNGRSP+PF2DTH-HFXOT-B 1 XOT- YNGXDT FFEGE X+r MORT- DE TRSP ) 

DE  TRSP=KJ*T2NC*AEROBCFKK*T2NC* ANOX IC 


NOSORT 

•THE  FOLLOWING  EQUATIONS  SCALE  THE  PRINTOUT  SO  THAT  THE  MINIMUM  IS  ZERO 

IF (TIME- 10*) 1 .2.2 
F ISH=ALLF 
FI SHNO=FN 
F I SH  0I  = FB 
V A SC  = VP 
PH YTO =PP 
ALGAE=EP 
I NVERT  = B I 
ZOOPL=ZP 
HERBF=HF 
PR  EOF 1 = PF 1 
PREDF2=PF2 

YOUNG=YNG 
DETR  TS=  DT 
PHOS 1 =PEP I 
PHCS2—PHYPO 
ALLPF1=ALLHF 
GO  TO  3 
F I SH=0. 

F ISHBI=0. 

F ISHNO=0. 

VASC=0. 

PHYT  0=0 . 

ALGAE=0 . 

I N VE  R T = 0 . 

ZOOPL=0. 

HERBF=0 . 

PREDF 1=0. 

PREOF2=0. 

YOUNG=0 . 

0ETRTS=0. 

PHOS 1=0. 

PHOS2=0 . 

ALLPF 1=0 
CONT INUE 


I 


- - r~ - — — Jl_  II  I "J 


- 


■ 


APPENDIX  B: 


EXPLANATION  OF  TERMS  USED  IN  THE  PROGRAM* 


Term  in  Equation  in  Figure  12 


Term  in  Computer  Program 


Phytoplankton  (Q^) 


Ku  • S 


(1  - R) 


13 


K1PP*PPL1TE 


hi. 


1 ♦ S • (1  - R)  • K, 


• f2Q2  • K. 


15 


PPMEFF 

i + PPLITE/PP50 


27 


fY  • (product  of  above  terms) 


CPPRL5*GPPPP 
( PPLEAC ) 


KUPP*T2PP*PP 
( PPRESP ) 


V KU3  ‘ V T5 


KHF2*PP*HF*FSHTMP 
( HFXPP) 

KZP2*PP*ZP*ZPTEMP 

(ZPXPP) 


“isA 


KYNG2*PP*YNG#FSH'IMF 
( YNGXPP ) 

k6pp*pp  (SINKNG) 


Explanation 


Availability  of  light 


Photosynthetic  rate  at  avail- 
able light  level 


Level  of  orthophosphate  in 
epilimnion 

Effect  of  temperature  on 
phot osynt hes i s 

Conversion  from  kilocalories 
to  carbon  for 
photosynthesis 

Phytoplankton  biomass 


Leaching  of  organic  carbon; 
rate  varies  with  rate  of 
photosynthesis 

Respiration  rate 


Consumption  by  herbivorous 
fish 

Consumption  by  zooplankton 


Consumption  by  young  primary 
predator  fish 

Settling  of  phytoplankton  to 
sediments 


Term  in  Equation  in  Figure  11 Term  in  Computer  Program Explanation 


Macrophytes  and  Epiphytic  Algae  (Q^)  and  Tubers  (Q^  ) 


R)  • 

-(K-Z'+K  Z- 

ve 

• Q^/Z) 

KJ  VP*VPLITE 

Availability  of  light 

VPMEFF 

Photosynthetic  rate  at  avail- 
able light  level 

i + ; 

LQa/t92.50 

1 ♦ VPLITE/VP50 

VP 

Macrophyte  biomass 

Tio 

T1VP 

Effect  of  temperature  on 
photosynthesis 

*27 

KCAL2C 

Conversion  from  kilocalories 
to  carbon  for 
photosynthesis 

fY  • 1 

S132  ‘ ^13 

BIRTH*KTU2*TU 

(SPROUT) 

Tuber  germination  in  spring 

(Continued) 

* A program  listing  for  the  Lake  Conway  model  is  given  in  Appendix  A. 

B1 


J ...  rru^r.un Explanation 


Macrophytes 

and  Epiphytic  Aigne  ((^  .uid  Tubers  (Q,  ,) 

(Continued) 

K24  ' S * *U 

K3VP*VP*T2VP 

(VPRET.P) 

Respiration  rate 

K25  ' % - « • K26  • s 

K6VP*VP  * K9VP*VP*BIE 
(SLOUGH) 

Normal  sloughing  rate  of 
leaves,  with  additional 
fall  pulse 

K131  • «2  ■ " 

DIE*KTUi»VP  (PRPCRH) 

Tuber  formation  in  fall 

V^VS 

ICBI2*BI*vrp*BrTEMP 
( BIXVP ) 

Consumption  by  benthic 
invertebrates 

®2  * K6!t " T1  ‘ ^6 

KHF3*VP*HF*FSHTMP 

(HFXVP) 

Consumption  by  herbivorous 
fish 

S ■ K8A  • T!  • «8 

KYNC  3*VP*YNG*FSHTMP 
( YNGXVP ) 

Consumption  by  young  primary 
predator  fish 

: .n  in  Equation  in  Figure  n 

Term  in  Computer  Program 

Explanation 

Epipelic  Algae  (Q^) 

• a-.,.  e'(K^  ‘ V 

• f3«2 

K1EP*EPLITE 

Availability  of  light 

A 

EPMEFF 

1 * «3/Ui3-50 

1 + EPLITE/EP50 

Photosynthetic  rate  at  avail- 
able light  level 

S.2 

PHYPO 

Level  of  orthophosphate  in 
hypolimnion 

T8 

T1EP 

Effect  of  temperature  on 
photosynthesis 

*27 

KCAL2C 

Conversion  from  kilocalories 
to  carbon  for 
photosynthesis 

s 

EP 

Biomass  of  epipelic  algae 

• (product  of  above  terms) 

CEPHLS'GPPEP 

(EPLEAC) 

Leaching  of  organic  carbon; 
rate  varies  with  rate  of 
photosynthesis 

K33-V«3 

K2EP*T2EP*EP( EPRESP ) 

Respiration  rate 

V “3 

K5EP»EP(EPM0RT) 

Death  rate 

KZF1*EP»ZP*ZPTEMP 

(ZPXEP) 

Consumption  by  zooplankton 

Q3  ' *62  ’ V % 

KHF1  *EP»HF*FSHTMP 
(HFXEP) 

Consumption  by  herbivorous 
fish 

wvs 

KPFl^EP'PFl'FSKTMP 

(PF1XEP) 

Consumption  by  adult  primary 
predator  fish 

®3  ‘ *82  ‘ T1  ‘ Q8 

KYNGl*r3>«YNG»FSHTMP 

(YNGXEP) 

Consumption  by  young  ; r «%  -. 
predator  fish 

Term  in  Equation  in  Figure  16 

Term  in  Computer  Program 

Zooplankton  (Q^) 

v ((IV  V * <ki,3- v * <*«.■  v> 

ZPASSM*ZPFOOD 

Food  assimilate:  ?.  . 
ton  (epipeli  al«a» 
plankton,  an : 

K1.5  ’ V *5 

KZP^*Zr*ZPTEMP 
( ZPRESP) 

Respiration  rst* 

(Continued) 

B2 
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Term  in  Equation  in  Figure  11 


Term  in  Computer  Program 


Explanation 


Macrophytes  and  Epiphytic  Algae  (Q^)  and  Tubers  (Q^) 

(Continued) 

*»*VTu 

K3VT*VP*T2VP 

(VFRESF) 

Respiration  rate 

*25  • - " • K26  • «2 

K6VP*VP  ♦ K9VP*VP*DIE 
(SLOUGH) 

Normal  sloughing  rate  of 
leaves,  with  additional 
fall  pulse 

*131  * V « 

DIE*JCTUi*VP  ( PRFORM) 

Tuber  formation  in  fall 

VSs'VS 

KBI2»BI*VP»BITEMP 

(BIXVP) 

Consumption  by  benthic 
invertebrates 

KHF3*VP*HF*FSHTMP 

(HPXVP) 

Consumption  by  herbivorous 
fish 

«2  • K8L  • T1  • «8 

KYNC3*VP*YNG*FSHTMP 
( YNGXVP ) 

Consumption  by  young  primary 
predator  fish 

Term  in  Equation  in  Figure  13 

Term  in  Computer  Program 

Explanation 

Epipelic  Algae  (^) 

-(K  2*K  • Q ) 

K31  ■ S • (1  - R)  • e 2 7 1 . f 

3*2 

K1EP*EPLITE 

Availability  of  light 

\ 

EPMEFF 

Photosynthetic  rate  at  avail- 
able light  level 

1 * «3/U!3-50 

1 ♦ EPLITE/EP50 

“l2 

PHYPO 

Level  of  orthophosphate  in 
hypolimnion 

T6 

T1EP 

Effect  of  temperature  on 
photosynthesis 

*27 

KCAL2C 

Conversion  from  kilocalories 
to  carbon  for 
photosynthesis 

s 

EP 

Biomass  of  epipelic  algae 

fY  • (product  of  above  terms) 

CEPRLC*GPPEP 
( EPLEAC ) 

Leaching  of  organic  carbon; 
rate  varies  with  rate  of 
photosynthesis 

K33  ‘ T8  ‘ ®3 

K2EP*T2EP*EP( EPRESP) 

Respiration  rate 

K3k*  S 

K5EP*EP(EPM0PT) 

Death  rate 

fCZPl*EP*ZP*ZPTEMP 

(ZPXEP) 

Consumption  by  zooplankton 

Va-V* 

KHF1 •EP*HF#FSHTMP 
(HFXEP) 

Consumption  by  herbivorous 
fish 

s-s-vs 

KPFlU*EP#PFl*FSHTMP 

(PF1XEP) 

Consumption  by  adult  primary 
predator  fish 

S’ ^ 

KYRG1*EP*YRG*FSHTMP 

(YHGXEP) 

Consumption  by  young  primary 
predator  fish 

Term  in  Equation  in  Figure  l6 

Term  in  Computer  Program 

Explanation 

Zooplankton  ( ) 

V VV<(Sa*  V * (*WV  * lKw 

■v> 

ZPASSM*ZPFOOD 

Food  assimilated  by  zooplank- 
ton (epipelic  algae,  phyto- 
plankton, and  macrophytes) 

Kl*5  * aU  ’ T5 

KZP5*ZP*ZPTEMP 
( Z PR ESP ) 

Respiration  rate 

(Continued) 


B2 


» 

Term  in  Equation  in  Fi 

Zooplankton  (Q^)  (Continued) 

VSp-S-Vt6 

KBI 1 *B  I *ZP*B  I TOtP^/ZPrEMT 
( BIXZP) 

Consumption  by  benthic 
invertebrates 

^ • K7U  • • T1  • Tk 

KPFl 3*ZP*PF1*FSHTMP2/ZPTEMP 
(PF1XZP) 

Consumption  by  adult  primary 
predator  fish 

^ • K67  • Hg  • T1  • TU 

KHF6»HF»ZP»FSHTMP2/ZPTEMP 
( HFXZP ) 

Consumption  by  herbivorous 
fish 

S ■ K87  • % ■ T1  • \ 

KYI«06»YNG*ZF*FSHTMP2/ZPTEMP 

(YNGXZP) 

Consumption  by  young  primary 
predator  fish 

Term  in  Equation  in  Figure  18 

Term  in  Computer  Program 

Explanation 

Benthic  Invertebrates  (Q^) 

•SrS'V'^-VV  * (k53-S»  ♦ 

K^)  biassm»bifood 

Food  assimilated  by  benthic 
invertebrates  ( zooplankton , 
macrophytes,  and  detritus) 

K5U  • S • T6 

kbii<»bi*bitemp(biresp) 

Respiration  rate 

VK66‘ VVT3 

KHF5*HF*BI*FSHTMP2 /BITEMP 
(HFXBI ) 

Consumption  by  herbivorous 
fish 

% ‘ *72  ' V V T3 

KPF11*BI«PF1*FSHTMP2/BITEMP 

(PF1XBI) 

Consumption  by  adult  primary 
predator  fish 

% • k86  • «B  • Tr  T3 

K7NG5*YNG«BI*FSHTMP2/BITEMP 
( YWGXBI ) 

Consumption  by  young  primary 
predator  fish 

S-SrVVV^ 

KFF2U*BI“PF2*(FSHTMP2/BITEMP) “COVER 
( PF2XBI ) 

Consumption  by  secondary  pre- 
dator fish 

Term  in  Equation  in  Figure  19 

Terra  in  Computer  Program 

Explanation 

Herbivorous  Fish  (Qg) 

K6rV  V((K62-S>  + (K63-V  + iK6U 

*K65  * (K66-W  * (K67  ‘ ' T3)) 

• Qg)  HFASSM*HFF00D 

Food  assimilated  by  herbivor- 
ous fish  (epipelic  algae, 
phytoplankton,  macrophytes, 
detritus,  benthic  inverte- 
brates, and  zooplankton) 

"68- VT1 

KHF7*HF*FSHTMP  ( HFRESP ) 

Respiration  rate 

*93  ’ V V V 'S 

KPF22*HF#PF2#FSHTMP*C0VER 

(HFPRED) 

Consumption  by  secondary 
predator  fish 

Term  in  Equation  in  Figure  20 

Term  in  Computer  Program 

Explanation 

Primary  Predator  Fish:  Adult  (Q^)  and  Young  (Qg) 


*n,VTi,<,lV% -V  + (kt3,q2> 

+ <*SU  * *Tk>  * (KT5,S,) 

Kaa ' 98  ‘ Ti 

V" 

VVVV'S 


P1ASSM*PF1FD 


KYNG7*YNG*FSHTMP 
( GROWTH ) 


Food  assimilated  by  adult 
primary  predator  fish  (ben- 
thic invertebrates,  macro- 
phytes, zooplankton,  and 
epipelic  algae) 

Growth  of  young  fish  into 
adult  group 


KPF15*PF1*BIRTH/DELT( BREED)  Birth  of  young  fish 


KPF21#PF2#PF1*FSH‘IMP*C0VEH 
( PF1PRD ) 


KPFl6*PFl»FSHTMP( PF1RSP) 


Consumption  of  adult  primary 
predator  fish  by  secondary 
predator  fish 

Respiration  rate 


(Continued) 

B3 


>lanation 


Term  In  Equation  In  Figure  20 Term  In  Computer  Program Exj> 


Primary  Predator  FI  ah: 

Adult  (Q^)  and  Young  (^g)  (Continued) 

K8l.  V V***' V ♦ (Kg,  • *j)  . K8k  • dg) 

* (K85- V * “BS'S'V  * (K87-  W> 

YNGASM*YNGFD 

Food  assimilated  by  young 
predator  fish  (epipelic 
algae  * phytoplankton,  de- 
tritus, benthic  inverte- 
brates, and  zooplankton) 

K89  ’ ®8  ‘ T1 

KYNG8*YNG*FnHTMP(  YMGRSf ) 

Respiration  rate 

V VS'V'S 

KPF? 3* PF2* YNG*FSHTMP*COVER 
(YNGPRD) 

Consumption  of  young  primary 
predator  fish  by  secondary 
predator  fish 

Term  in  Equation  in  Figure  21 

Term  in  Computer  Program 

Explanation 

Secondary  Predator  Fish  (Q^) 

lsrvv((W  * ‘•VW'V 

* <«*’ VV*V  * (K95  • V ti  ‘ ’V  fqa)) 

PF2  ASM*  PF2  FD 

Food  assimilated  by  secondary 
predator  fish  (young  and 
adult  primary  predator  fish, 
herbivorous  fish,  and  ben- 
thic invertebrates) 

**•  VTi 

KPF23*PF?*FSHTMP 
( PF2RSP ) 

Respiration  rate 

Vr'V’i 

KPF?6*PF2*FSHTMP 

(PF2DTH) 

Death  rate 

Term  in  Equation  in  Figure  22 

Term  in  Computer  Program 

Explanation 

Detritus  (Q^q) 

■WV" 

KK*T2NC*AN0XIC 

Release  of  orthophosphate 
from  sediments  when 
hypolianlon  is  anaerobic 

KJ  *T2NC • AERO  BC 

Release  of  orthophosphate 
from  sediments  when 
hypolimnion  is  aerobic 

Term  in  Computer  Program 

Explanation 

Orthophosphate  in  Eplllmnlon  (Q-.) 

n’S 

0P1HFT 

Orthophosphate  entering  as 
rainfall  and  runoff 

K 

P 

KC2P 

Carbon  to  phosphorus  ratio  in 
organic  matter 

(«X2-«U>-" 

TRNOVR* ( PHYPO-PEPI ) 

Mixing  between  epilimnion  and 
hypolimnion  during 
isothermal  periods 

Term  in  Computer  Program 

Explanation 

White  Amur: 

Biomasa  (Q^)  And  Number  (Q^) 

Kllil  ‘ *1112  ‘ Klk3  ’ Sk  ‘ S ’ T1 

K1FB*K2FB*K3FB*FB*VP*FSHTMP 

(PASSIM) 

Assimilation  rate  of 
macrophytes 

Sm’Sk"! 

K5FB#rB»FSHTMP(  FRESF ! 

Respiration  rate 

*151  ' *15  ’ S>i/915 


t 


I 


1 


IQFH»FII*FIi/FII(nciRT) 


Mortality  rate 


In  accordance  with  letter  from  DAEN-RDC,  DAEN-ASI  dated 
22  July  1977,  Subject:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a facsimile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 
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